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ABSTRACT 

Over  the  past  two  decades,  materials  engineers  have  developed  piezoelectric-polymer 
composite  materials  that  enable  effective  electromechanical  properties  to  be  tailored  for 
a  specific  application.  These  materials  are  made  by  combining  conventional  piezoelectric 
ceramics  with  piezoelectrically  passive  polymers  in  a  variety  of  geometrical 
configurations.  As  with  any  composite  material,  the  properties  and  behavior  of 
piezocomposites  are  highly  dependent  on  the  properties  of  the  constituent  materials  and 
the  local  arrangement  of  the  different  phases.  In  particular,  the  ceramic-polymer 
interface  plays  an  important  role  in  determining  the  electromechanical  coupling  in  the 
piezocomposite.  In  this  dissertation,  the  electromechanical  behavior  of  1-3 
piezocomposites  is  investigated  from  both  theoretical  and  experimental  stand  points. 

Theoretical  investigations  centered  on  the  development  of  a  micromechanics  model 
for  predicting  the  local  fields  and  effective  behavior  in  piezocomposites  with  1-3 
connectivity.  Since  the  presence  of  a  thin  interlayer  or  polymer  coating  around  the 
ceramic  rods  can  influence  the  local  interaction  between  the  piezoceramic  and  polymer 
matrix  and  change  the  overall  performance  of  the  composite,  a  finite  composite  cylinder 
model  was  developed  to  incorporate  an  interlayer  with  varying  properties.  Experimental 
studies  focused  on  probing  the  surface  displacements  of  1-3  piezocomposites  using  a 
scanning  heterodyne  laser  interferometer.  Static  surface  displacements  of  1-3  PZT  rod- 
epoxy  samples  with  different  interphase  regions  were  measured  and  correlated  with  the 


effective  low-frequency  performance  of  the  composite.  Several  types  of  interphase 
region  were  considered.  Coatings  with  elastic  moduli  lower  than  that  of  the  epoxy  matrix 
were  applied  to  the  rods.  The  influence  of  a  silane  coupling  agent  was  also  investigated. 
Experimental  displacement  profiles  were  compared  with  micromechanical  predictions 
using  the  finite  composite  cylinder  model.  The  results  demonstrate  that  the  presence  of 
an  interphase  between  the  piezoceramic  and  the  polymer  matrix  influences  the  local 
deformations  and  changes  the  overall  performance  of  the  composite.  Thus,  the  interphase 
plays  an  important  role  in  determining  the  electromechanical  coupling  in  the 
piezocomposite. 

The  study  of  electromechanical  coupling  in  piezocomposites  was  further  developed 
by  investigating  the  hydrostatic  performance  of  1-3  piezocomposites,  an  important  issue 
in  design  of  piezocomposites  for  low-frequency  applications.  Emphasis  was  placed  on 
determining  the  stress  transferred  between  the  passive  matrix  and  the  active  piezoceramic 
rods  and  using  this  data  to  indicate  the  level  of  electromechanical  coupling.  The  stress 
field  in  the  piezoelectric  ceramic  under  hydrostatic  loading  was  predicted  using  the 
analytical  micromechanical  model  developed  and  a  finite-element  model  as  well. 
Optimal  electromechanical  coupling  was  achieved  when  a  certain  favorable  stress  field 
was  induced  in  the  piezoceramic.  The  influence  of  such  design  parameters  as  the  matrix 
stiffness,  the  interphase  stiffness,  the  interphase  thickness,  the  Poisson’s  ratio  of  the 
polymers,  and  piezoceramic  rod  aspect  ratio  on  the  hydrostatic  performance  of  1-3 
piezocomposites  was  also  investigated.  Although  the  current  work  is  focused  on  the 
electromechanical  behavior  of  1-3  piezocomposites  at  low-frequency,  the  research  results 
and  conceptual  understanding  obtained  have  importance  for  optimizing  the  design  of 
piezocomposites  in  other  applications  as  well. 
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1.  INTRODUCTION 


1.1  Piezoelectricity 

Piezoelectricity,  according  to  the  original  definition  of  the  phenomenon  discovered  by 
Jacques  and  Pierre  Curie  brothers  in  1880  (Jaffe,  Cook  and  Jaffe,  1971),*  is  the  ability  of 
certain  crystalhne  materials  to  develop  an  electric  charge  which  is  proportional  to  a  mechanical 
stress.  Soon  after  this  discovery,  it  was  also  realized  that  the  converse  effect  of  a  geometric 
strain  developing  upon  the  apphcation  of  a  voltage,  was  also  inherent  in  these  materials.  Thus, 
piezoelectricity  is  an  interaction  between  electrical  and  mechanical  systems.  The  direct 
piezoelectric  effect,  the  development  of  an  electric  charge  upon  the  application  of  a  mechanical 
stress,  is  described  as  (Nye,  1957) 

Pi=d,Jj,  (1.1) 

where  is  a  component  of  the  electric  polarization  (charge  per  unit  area),  are  the 
components  of  piezoelectric  coupling  coefficient,  and  T are  components  of  the  apphed 

mechanical  stress.  The  converse  effect,  the  development  of  a  mechanical  strain  upon  the 
apphcation  of  an  electric  field  to  the  piezoelectric,  is  described  by 

=  (1.2) 

where  S^j  is  the  strain  produced  and  E^.  is  the  applied  electric  field.  In  both  cases,  the 
piezoelectric  coefficients  d.j^.  are  numericaUy  identical. 

The  piezoelectric  effect  is  strongly  linked  to  the  crystal  symmetry.  Piezoelectricity  is  limited 
to  20  of  the  32  crystal  classes,  known  as  point  groups,  for  all  crystalhne  materials.  The  crystals 


References  are  arranged  alphabetically  by  author  in  the  Bibliography,  beginning  on  page  102. 
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which  exhibit  piezoelectricity  have  one  common  feature:  the  absence  of  a  center  of  symmetry 
within  the  crystal.  This  absence  of  symmetry  leads  to  polarity,  the  one-way  direction  of  the 
charge  vector.  Most  of  the  important  piezoelectric  materials  are  also  ferroelectric,  which  is 
characterized  by  a  transformation  to  high  symmetry,  non- piezoelectric  phase  at  higher 
temperatures.  The  transformation  temperature  is  known  as  the  Curie  temperature.  For  any 
piezoelectric  materials,  the  Curie  temperature  is  the  absolute  maximum  use  temperature. 

PolycrystaUine  materials  in  which  the  crystal  axes  of  the  grains  are  randomly  oriented  exhibit 
zero  net  piezoelectric  effect.  Piezoelectric  polycrystaUine  ceramics  were  discovered  in  the 
1940s,  soon  followed  by  the  development  of  the  poling  process  in  which  the  randomly  oriented 
crystal  axes  are  suitably  aligned  by  the  apphcation  of  a  strong  electric  field  at  elevated 
temperatures.  After  that  discovery,  researchers  continued  to  develop  better  and  more  stable 
materials.  Lead  zirconate  titanate  (PZT)  was  first  introduced  in  1954  and  has  become  the  most 
widely  used  piezoceramic.  Today  piezoelectric  ceramics,  or  piezoceramics,  are  used  more 
widely  than  piezoelectric  crystals. 

Generally,  mechanical  stress  T-j  and  electric  field  E-  are  coupled  within  a  piezoelectric  sohd 
medium.  Constitutive  relations  for  a  piezoelectric  material  can  be  written  as  (IEEE,  1978) 

^ij  ~  ^ijkl  ^kl  ^ijk^k  (1-3) 

A  =  +  (1-4) 

where  D,  is  the  dielectric  displacement  given  by  the  sum  -I-P,.  Here  8^  are  dielectric 
constants,  are  elastic  comphance  components,  and  eg  is  the  dielectric  constant  of  free 
space.  Superscripts  T  and  E  have  been  added  to  8^  and  to  show  that  these  constants 
describe  dielectric  and  elastic  properties  measured  under  conditions  of  constant  stress  and 
constant  electric  field,  respectively. 
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The  constitutive  relations  can  be  written  in  alternate  forms,  such  as 

~  ^ijkl  ^kl  ~  ^ijk ^k  (1-5) 

D,=e,^,S^,  +  4E^.  (1.6) 

where  the  piezoelectric  constants  e^j,^  are  related  to  the  piezoelectric  constants  and  the 
elastic  stiffness  constants  by  the  expression  An  abbreviated  subscript 

notation  can  be  introduced  to  reduce  the  three  subscripts  of  the  piezoelectric  constants  and  the 
four  subscripts  of  the  elastic  constants  into  two  subscripts.  Relations  between  the  piezoelectric 
constants  and  elastic  constants  with  full  subscripts  and  those  with  abbreviated  subscripts  are 
outlined  in  the  IEEE  Standard  on  Piezoelectricity  (IEEE,  1978). 

In  Eqs.  (1.3)  and  (1.4)  the  strain  and  electric  field  are  derivable  from  the  mechanical 
displacement  and  electric  potential  as 


(1.7) 

(1.8) 

To  complete  the  formulation  of  the  stationary  hnear  theory  of  piezoelectricity,  one  must 
supplement  Eqs.  (1.3)  to  (1.6)  with  the  equations  of  elastic  equihbiium  and  the  electrostatic 
equation,  which  in  the  absence  of  body  forces  and  free  electric  charge  are 

=  Q  (1.9) 

D,,  =  0.  (1.10) 


3 


Piezoelectric  materials  are  widely  used  in  hydro-  and  electro-acoustics,  electro- optics, 
communications,  and  measurement  techniques  (Dokmeci,  1980).  The  research  on  piezoelectric 
materials  and  piezoelectricity  has  resulted  in  a  huge  volume  of  literature.  Two  references  that 
provide  an  overview  on  this  topic  are  included  in  the  Bibliography  (IEEE,  1978;  Ikeda,  1990). 

1.2  Piezoelectric  Composites 

1.2.1  Background 

Composite  piezoelectric  materials  have  been  developed  over  the  past  the  two  decades 
which  enable  the  effective  properties  to  be  tailored  for  a  specific  application.  These  materials 
are  made  by  combining  conventional  piezoelectric  ceramics  with  piezoelectrically  passive 
polymers  in  a  variety  of  geometrical  configurations.  Piezocomposites  were  originally  developed 
for  underwater  hydrophone  apphcations  in  the  low-frequency  range,  but  have  also  been 
extended  to  such  other  apphcations  as  ultrasonic  transducers  for  acoustic  imaging  and  medical 
apphcations  (Gumraja,  Safari  et  al.,  1988)  .  Recently,  piezoelectric  ceramics  are  being 
investigated  for  use  as  sensor  or  actuator  elements  in  smart  stmctures. 

Unlike  stmctural  composite  materials  which  are  formulated  to  optimize  mechanical 
properties,  piezocomposites  are  designed  to  maximize  their  coupled  field  behavior.  The  key 
consideration  in  designing  a  piezocomposite,  whether  as  a  sensor  or  as  an  actuator,  is  to  ensure 
the  maximum  efficiency  in  electromechanical  energy  conversion  (Smith,  1992).  As  with  any 
composite  material,  the  properties  and  behavior  of  piezocomposites  are  highly  dependent  on  the 
properties  of  the  constituent  materials  and  the  local  arrangement  of  the  different  phases.  In 
particular,  the  ceramic -polymer  interface  plays  an  important  role  in  determining  the 
electromechanical  couphng  in  the  piezocomposite  (Li  and  Sottos,  1995). 

The  first  attempt  at  making  piezoelectric  composites  occurred  in  the  early  1970s  when 
researchers  tried  to  develop  flexible  piezoelectrics  (Pauer,  1973).  In  the  late  1970s  a 
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comprehensive  research  program  on  piezoceramic -polymer  composites  was  carried  out  at 
Pennsylvania  State  University.  Piezocomposites  have  since  become  an  important  addition  to  the 
range  of  composite  materials,  and  numerous  papers  have  been  written  regarding  manufacturing, 
modeling,  vibration  modes,  and  various  apphcations.  Piezocomposites  have  advantages  over 
single-phase  piezoelectrics  for  certain  apphcations.  Gumraja  et  al.  summarized  the  pioneering 
work  on  piezocomposites  at  Pennsylvania  State  University  in  a  section  of  the  book  Electronic 
Ceramics  pubhshed  in  1988  (Gururaja,  Safari  et  al.,  1988).  A  review  article  of  ah  work  prior 
to  1989  was  presented  by  Smith  at  the  1989  TERR  Ultrasonic  Symposium  (Smith,  1989). 

A  variety  of  piezocomposites  can  be  made  by  combining  piezoelectric  ceramics  with 
passive  polymers.  These  composites  are  commonly  classified  according  to  their  connectivity 
(Newnham,  1986).  Connectivity  in  a  composite  refers  to  the  way  in  which  the  different  phases 
are  arranged,  and  is  designated  by  such  numbers  as  1-3,  3-2,  etc.  The  numbers  represent  the 
number  of  dimensions  in  which  each  phase  is  self- connected.  Some  piezocomposites  with 
different  connectivities  are  shown  in  Fig.  1.1. 
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(d) 


Fig.  1.1.  Piezocomposites  with  different  connectivities. 

(a)  A  1-3  piezocomposite  with  piezoceramic  rods  embedded  in  a 
continuous  polymer  matrix. 

(b)  A  2-2  piezocomposite  found  in  multilayer  piezoceramic  actuators. 

(c)  A  piezocomposite  with  0-3  connectivity. 

(d)  The  moonie  2-2  piezocomposite,  an  advanced  piezoelectric  sensor 
and  actuator. 
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The  literature  regarding  piezocomposites  can  be  roughly  divided  into  five  basic  categories: 
(1)  micromechanical  analysis  of  low-frequency  effective  elastic,  dielectric,  and  piezoelectric 
behavior  and  use  of  low-frequency  parameters  in  hydrophone  and  one -dimensional  transducer 
models;  (2)  prediction  of  the  various  resonance  modes  associated  with  the  composites  and  their 
coupling  to  one  another  (including  experimental  and  finite- element  analysis  studies);  (3) 
presentation  of  techniques  for  composite  manufacture;  (4)  reports  of  composite  transducer 
performance  results  and  (5)  reports  on  smart- material  and  smart  stmcture  systems  apphcations. 
Papers  written  on  the  first  topic  are  most  critical  for  the  current  investigation  and  are 
summarized  in  the  following  section. 

1.2.2  Micromechanical  Models 

Most  of  the  existing  micromechanics  works  concerning  the  electro- mechanical  coupling 
phenomena  in  piezocomposites  involve  the  derivation  of  relations  for  the  effective  elastic 
constants  and  effective  piezoelectric  constants.  Classical  models,  such  as  the  Voigt  uniform- 
strain  model  and  the  Reuss  uniform- stress  model,  have  been  utihzed  (Lees  and  Davidson, 
1977).  More  refined  models  include  the  parallel- serial  model  (Haun  and  Newnham,  1986)  and 
the  model  by  Smith  (1993).  These  simple  models  are  aU  based  on  assumption  that  either  strain 
components  or  the  corresponding  stress  components  are  equal  in  the  two  phases.  In  addition, 
Jensen  (1991)  has  used  a  concentric  cylinder  model  to  characterize  1-3  piezocomposites.  This 
model,  which  estimates  the  six  technically  most  important  constants  of  the  nine  constants  in  the 
e-set,  consists  of  a  piezoelectric  rod  and  a  concentric  elastic  tube  subjected  to  axisymmetric 
load.  Because  the  displacement  assumption  introduced  in  Jensen’s  study  is  equivalent  to  that  of 
generahzed  plane  strain,  the  concentric  tube  model  also  leads  to  homogeneous  strain  along  the 
axial  direction  of  the  cy Under.  Cao,  Zhang  and  Cross  (1992)  presented  a  theoretical  study  on 
the  static  performance  of  1-3  piezocomposites.  Inhomogeneous  displacement  profiles  were 
derived  for  a  single-rod  composite  and  a  single-tube  1-3  ceramic -polymer  composite  under 
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uniaxial  or  hydrostatic  stress.  The  effect  of  lateral  stress  is  not  included  in  the  study  and  the 
analysis  is  essentially  one- dimensional. 

Recent  developments  of  micromechanics  of  piezocomposites  also  include  the  work  of  Dunn 
and  Taya  (1993),  who  estimated  the  effective  properties  of  a  two-phase  piezocomposite  using 
dilute,  self- consistent,  Moii-Tanaka  and  differential  micromechanical  models.  Progress  in 
micromechanical  analysis  of  piezocomposites  has  also  been  made  by  Benveniste  (1993a, 
1993b)  and  Benveniste  and  Dvorak  (1992).  Benveniste  and  Dvorak  (1992)  extended  the 
uniform- fields  concept  for  the  uncoupled  mechanical  case  to  the  piezoelectric  behavior  in  binary 
piezocomposites  and  derived  the  universal  relations  for  the  effective  piezoelectric  constants  of 
such  composites.  More  recently,  Benveniste  (1994a)  combined  the  composite  cylinder 
assemblage  model  with  these  universal  relations  for  the  effective  piezoelectric  constants  and 
derived  some  exact  expressions  for  the  effective  constants  that  characterize  two-phase  fibrous 
piezocomposites.  A  summary  and  discussion  of  his  recent  contribution  in  the  micromechanics  of 
piezocomposites  was  given  by  Benveniste  recently  (Benveniste,  1994b). 

1.2.3  Experimental  Characterization 

Experimental  characterization  of  piezocomposites  can  be  roughly  divided  into  two 
categories.  The  first  is  associated  with  measurements  of  dielectric  and  piezoelectric  constants 
by  using  either  the  direct  piezoelectric  effect  or  the  resonance  method  outlined  in  the  IEEE 
Standard  on  Piezoelectricity  (IEEE,  1978).  The  second  involves  measuring  surface  velocity  and 
displacement  at  certain  frequencies  by  laser  probing  methods. 

Klicker  (1980)  used  a  Berlincourt  piezo  J33  meter  and  measured  the  J33  value  at  different 

locations  on  the  surface  of  a  1-3  piezocomposite.  He  then  chose  the  average  value  as  the 
measured  effective  piezoelectric  constant  of  the  composite.  This  method  was  also  used 
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recently  by  Kim,  Rittenmyer  and  Kahn  (1993).  Klicker  (1980)  also  used  a  static  technique 
based  on  the  direct  effect  to  measure  the  effective  hydrostatic  piezoelectric  constant  . 

Measuring  surface  displacement  at  resonance  frequencies  using  a  laser  heterodyne 
interferometer  was  reported  by  Wang  and  Auld  (1985)  and  Gumraja,  Schulze  et  al.  (1985). 
Rittenmyer  and  Dubbelday  (1992)  reported  a  direct  measurement  of  the  temperature- 
dependent  piezoelectric  coefficients  of  composite  materials  by  laser- Doppler  vibrometry. 
Zhang  et  al.  (1994)  measured  the  surface  displacement  profile  of  a  2-2  piezocomposite  at  low 
frequency  (200  Hz)  by  a  double -beam  laser  dilatometer.  The  low-frequency  performance  of  a 
single  rod  piezocomposite  has  been  investigated  by  Sottos  et  al.  (1993)  by  a  laser  heterodyne 
interferometer.  However,  characterizing  the  static  deformation  in  piezocomposites  has  not  been 
presented  in  the  literature. 

1.3  Overview  of  the  Dissertation 

This  dissertation  investigates  the  behavior  of  electromechanical  couphng  in  piezocomposites. 
The  work  focuses  on  several  topics  of  practical  importance  that  have  not  been  adequately 
addressed  in  the  existing  literature.  These  topics  include  the  load  transfer  between  the 
piezoceramic  and  the  polymer  matrix;  the  micromechanics  of  piezocomposites,  incorporating 
any  existing  interphase  between  the  ceramic  and  the  matrix;  the  optimization  of 
electromechanical  coupling  by  varying  interphase  properties;  and  the  experimental 
characterization  of  the  performance  of  piezocomposites  for  low  frequency  apphcations.  Special 
attention  has  been  given  to  the  important  issue  of  tailoring  the  ceramic -polymer 
interface/interphase  for  optimal  performance.  The  study  has  been  carried  out  from  both 
theoretical  and  experimental  standpoints.  The  balance  was  achieved  by  conducting  innovative 
laser  interferometric  measurements  and  developing  analytical  and  numerical  (finite- element) 
models  to  study  the  performance  of  the  materials.  Although  the  work  is  focused  on  the 
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performance  of  1-3  piezocomposites  at  low  frequency,  the  research  results  have  importance  for 
other  applications  as  well. 

In  Chapter  2,  a  concentric  composite  cylinder  model  is  developed  to  study  the  local  fields 
and  effective  behavior  in  piezoceramic  rod- polymer  composites  with  1-3  connectivity.  This 
composite  model  is  similar  to  the  well  known  “composite  cyhnder  assemblage”  model,  which  is 
used  commonly  for  continuous  fiber  composites.  This  composite  model  is  characterized  by  its 
finite  length,  chosen  to  equal  to  that  of  the  embedded  piezoceramic  rods.  The  displacement, 
stress,  and  electric  fields  in  the  composite  model  are  found  using  the  displacement  and  electrical 
potentials.  The  solution  is  compared  with  solutions  of  simple  plane  strain- models  and  other 
micromechanical  models.  The  effective  piezoelectric  constants  ^33  and  for  the  composite 

are  then  calculated  using  the  derived  solution. 

Chapter  3  focuses  on  experimental  characterization  of  local  deformation  in  1-3 
piezocomposites.  A  laser  heterodyne  micro- interferometer  is  used  to  study  the  static 
performance  of  1-3  piezocomposites.  Surface  displacement  profiles  of  composite  samples  are 
measured  around  the  piezoceramic  rods.  The  displacement  profiles  are  also  compared  with  the 
analytical  results. 

Electromechanical  couphng  is  discussed  further  in  Chapter  4.  The  focus  is  on  optimizing 
hydrostatic  performance  of  1-3  piezocomposites,  an  important  issue  in  design  of 
piezocomposites  for  low-frequency  apphcations.  Emphasis  was  placed  on  determining  the 
stress  transferred  between  the  passive  matrix  and  the  active  piezoceramic  rods  and  using  this 
data  to  indicate  the  level  of  electromechanical  couphng.  The  stress  field  in  the  piezoelectric 
ceramic  is  studied  using  the  solutions  developed  in  Chapter  2.  Optimal  electromechanical 
coupling  can  be  achieved  when  a  certain  favorable  stress  field  is  induced  in  the  piezoceramic. 
The  influence  of  such  design  parameters  as  the  matrix  stiffness,  the  interphase  .stiffness,  the 
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interphase  thickness,  the  Poisson’s  ratio  of  the  polymer  and  piezoceramic  rod  aspect  ratio,  on 
the  hydrostatic  performance  of  1-3  piezocomposites  is  investigated. 

Finally,  Chapter  5  summarizes  the  important  findings  of  the  thesis.  Concluding  remarks  and 
recommendations  for  future  work  are  presented. 


11 


2.  MICROMECHANICS  OF  PIEZOCOMPOSITES 


2.1  Introduction 

Consider  a  1-3  piezocomposite  with  an  array  of  aligned  piezoelectric  ceramic  rods 
embedded  in  a  passive  polymer  matrix,  as  shown  in  Fig.  2.1.  The  piezoelectric  ceramic  rods 
are  a  poled  polycrystahine  ceramic  which  is  isotropic  in  a  plane  normal  to  the  pohng  direction 
and  effectively  exhibits  transversely  isotropic  elastic  behavior.  Piezoelectric  ceramics  of  this 
family,  such  as  barium  titanate,  lead  zirconate  titanate  (PZT),  lead  metaniobate,  and  modified 
lead  titanate,  have  been  developed  for  a  broad  spectrum  of  apphcations.  The  convention  is  to 
define  the  pohng  direction  as  the  direction  and  the  x^-  X2  plane  as  the  isotropic  plane. 
Constitutive  relations  for  piezoceramics  poled  in  the  x,  direction  can  then  be  written  in  the 
matrix  form  (IEEE,  1978) 
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Ceramic  rod  Polymer  matrix 


Fig.  2.1.  Schematic  of  a  1-3  piezocomposite. 
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The  concentric  finite  composite  cylinder  model  shown  in  Fig.  2.2  is  chosen  as  the 
representative  volume  element  (RVE)  for  the  1-3  piezocomposite.  The  length  of  the  composite 
cylinder  model  is  chosen  as  the  length  of  the  piezoceramic  rods  embedded  in  the  polymer 
matrix.  If  the  length  of  the  rods  approaches  infinity  i.e.,  if  the  rods  are  continuous,  this 
composite  cylinder  model  wiU  approach  the  well  known  “composite  cylinder  assemblage”  of 
Hashin  and  Rosen  (1964)  and  Hashin  (1979). 

Because  the  actual  stress  transfer  is  through  the  interface  between  the  rods  and  the  matrix, 
the  incorporation  of  any  existing  interphase  region  into  the  micromechanical  analysis  of  1-3 
piezocomposites  may  be  important  for  understanding  composite  behavior.  A  thin  interlayer  is 
introduced  between  the  piezoceramic  rods  and  the  matrix  in  the  micromechanical  model.  It  is 
proposed  that  the  presence  rf  this  thin  interlayer  can  influence  the  effective  behavior  of  the 
composite  and  change  the  sensitivity.  Both  the  interlayer  and  the  matrix  are  considered  as 
isotropic  and  nonpiezoelectric. 

2.2  Formulation  of  the  Micromechanical  Analysis 

As  shown  in  the  schematic  diagram  of  the  three-phase  composite  cylinder  in  Fig.  2.2,  the 
region  r  <r^  (ceramic  rod),  the  region  r^<r<r^  (interlayer)  and  the  region  <r  <r^ 

(embedding  polymer  matrix)  are  denoted  by  superscripts  (1),  (2)  and  (3),  respectively.  The 
composite  aggregate  is  subjected  to  homogeneous  boundary  conditions  in  the  form 


Pi(^)  =  T°n^  (2.3) 

(^(f2)=-E;’x,  (2.4) 

where  p.  denotes  the  traction  vector,  (|)  is  the  electric  potential,  T°  is  a  constant  stress  tensor, 
Ff  is  a  constant  electric  field  vector,  n.  is  the  normal  of,  ^2,  the  outer  surface  of 

‘  J 


14 


Fig.  2.2.  Schematic  of  the  three-phase  composite  cyhnder  model. 
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the  composite  cylinder,  and  x,  denotes  the  components  of  the  Cartesian  system.  Because  of 
the  symmetry  of  the  geometry  and  the  boundary  conditions,  the  deformation  is  also 
axisymmetric. 


In  cylindrical  coordinates,  the  axially  symmetric  state  of  deformation  is  characterized  by 


=  u 

\r,z). 

(2.5) 

Mg  ^  =  0 

(2.6) 

ul  -  w 

'\r,z) 

(/  =  1,2,3). 

(2.7) 

For  a  transversely  isotropic  piezoelectric  rod  he  constitutive  relations  can  be  written  in 
cylindrical  coordinates  as 


'T’(l) 

1 


7^(1) 
^  06 


'T’(l) 
^  zz 

T^(l) 
^  rz 


^^12^00  ^^n^zz  ^3\ 


(1)77(1) 


-  ^(Dcd)  .  ^(l)c(l)  .  ^(1)C(1)_  ^(1)/7(1) 
“^12^rr  ^^11^06  ^^n^zz  ^3\^z 


-  ^a)c(i)  .  ^(i)c(i)  .  ^(i)c(i)_  ^(i)/7(i) 
~^n^rr  ^^13^06  ^33  ^z 


'15  ' 


y^(l)  _  y^(l)  _  ^ 

''re  ~  ~  ^ 


(2.8) 

(2.9) 

(2.10) 

(2.11) 

(2.12) 


(l)c(l) 


,(1)^(1) 


De®  =0 


n(l)  _  .(l)c(l)  ,  .(Dcd)  .  .(l)c(l)  ,  p(l);7(l) 

^Z  ~^3\'^rr  ^^31^6  ^^33^zz  ^^33^z  * 


(2.13) 

(2.14) 

(2.15) 


The  constitutive  relations  for  the  isotropic,  nonpiezoelectric  interlayer  and  matrix  have  the 
form 


(0c(‘) 

ZZ 


7^(0  _  ^(0c(0  ,  ^(0c(0  ,  ^(U  c 

1  rr  “  ^11  ^rr  ^  ^12^06  ^  ^12^; 

7^(0  _  ^(0c(()  _|_  ^(0c(0  _|_^(0c(0 
^  06  “  ^12^rr  ^11  ^60  ^12^; 

7^(0  _  ^(0c(0  ,  ^(0c(0  ,  ^(U  c 
^  zz  ~^\2^rr  ^^12^60^^11^; 


zz 
(0c(0 


(2.16) 

(2.17) 

(2.18) 
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a  =  2,3) . 


(2.19) 


T 


(0  - . 


•(q«-c«)5 


(0 

rz 


Substitution  of  the  constitutive  relations  and  the  strain- displacement  relations  into  the 
equihbiium  equations  and  the  electrostatic  equation  gives  the  following  governing  equations  in 
terms  of  displacements  and  electric  field  (/=1,2,3): 


JO 

'11 


f3  2,.(o 


.in 


-  e 


dr 

”  3r 


r  dr 


+k:+C) 


w 


(,) 

^  drdz  az^ 


-  e 


0  dE 


(0 


dz 


0 


(2.20) 


JO 

'44 


f  a  ^  1  aw«  ^ 


-I-  c 


dr^ 

av''^ 


ar 


+  (Ci3  +  C44) 


drdz 


■-I-- 


1  du 


(0  A 


az 


-I- 


ji, 


3z" 


'33 


az 


'15 


JO  J 


dr 


=  0 


(2.21) 


JO 

'15 


'a  ^  1  aw^'^ ' 

dr^  r  dr 
V  J 


-I-  e 


a 


dz 


l^Cl  +  p(0^^  ,  p(0 

“•"^33  “'"^11 


dz 


^  dEl‘^ 

dr 


^a^.jo 


drdz 


+  ■ 


1  du 


(0  J 


dz 


-I- 


jo  J 


=  0 


(2.22) 


For  the  ceramic  rod  (region  1),  the  governing  equations  for  mechanical  deformation  and 
electric  field  are  coupled.  The  equations  are  solved  by  introducing  two  displacement  potential 
functions  and  one  electrical  potential  function  such  that 

(1),  ,  at/(r) 

u  (r,z)  =  — —  cospz  (2.23) 

dr 

w^^\r,  z)  = -W {r)sm  \\z  (2.24) 
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(2.25) 


\r,z)  =  — — sm  |iz 
or 


E'f  V,  z)  =  |l£'(r)cos  |iz 


(2.26) 


where  U,  W  and  E  are  functions  of  r  only. 


In  order  to  satisfy  the  equations  of  equihbiium  and  the  electrostatic  equation,  the  functions 
U,  W  and  E  must  satisfy  the  equations 


(2.27) 

(2.28) 

(2.29) 

where  V j  =  — ^  . 

or  rar 


(4Xi  -  w  -  (c  ® )^w  -  )ia£ = 0 

(ciV  +  4'  wp + (PIP  -  PIP ) w + ) ^  =  0 

(63^;)  +  els'  WP  +  (elPl  -  4V"  )W  -  (el?vf  -  8''3V4£=  0 


Equations  (2.27),  (2.28),  and  (2.29)  can  be  further  reduced  into  a  higher  order  ordinary 
differential  equation 

(EP  +  hP  +  hP  +  h) {U,  W,E}=Q  (2.30) 

where  Pj ,  P2,  E,  and  are  dependent  on  the  material  properties  of  the  piezoelectric  rod 
and  are  given  in  the  Appendix. 

The  solutions  of  Eq.  (2.30)  for  a  solid  rod  are 

U  =  AEE(a^r)  +  A''Vo(a2r)  +  A^/,,  (a3r)  (2.3 1) 

W  =  A/'>P,/o(  a,r)  +  pph(^2r)  +  'P  34  (^^r)  (2.32) 
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E=  +  +  (2.33) 

where  a,  are  the  roots  of  the  equation 

Pittf  +P^a^^+P^a]+P^=0  (2.34) 

and  the  constants  (3,  and  y ,  are  given  by 

(ciV  +c®)l^P/  +(4?  =(Cn’ocf-|i"c'4^)  (2.35) 

(ciV« f  )P,.  +  (ej'^af  -  (t'e^s^) y,.  =  -(c^V  +  c®  )fia?  •  (2.36) 

The  structure  of  the  solutions  of  Eqs.  (2.31),  (2.32),  and  (2.33)  depends  on  the  nature  of 
the  value  of  af ,  and  hence  depends  on  the  values  of  the  material  properties.  For  a  PZT-5H 
rod  (Morgan  Matroc  Inc.,  1989),  a\  andaj  are  complex  conjugates.  Because  the 
displacements  and  stresses  are  real,  the  functions  U,  W  and  E  are  chosen  in  the  form 

U  =  A"4(air)  +  Re[/o(a2r)]+  Im[/o(a2r)]  (2.37) 

W  =  A<')p,/„(a,r)  +  A«  Re[p2/o(a20]+  Af >  Im[p2/o(a2^)]  (2-38) 

E  =  A«y/„(a^r)  +  A«  Re[y2/o(a2r)]  + A<^^  Im[y2/o(a2r)]  (2.39) 

where  4  is  the  zero- order  modified  Bessel  function  of  the  first  kind. 

The  complementary  displacements,  stresses  and  dielectric  displacements  for  the  ceramic 
rod  (/=1)  are  then  obtained  as 

uT  =  £cos(p,„z)[a?/2i?(?-)  +  A^^h2Xr)  +  (2.40) 
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sin(  |j,„ z)Wjh^^J(r)  +  (r)  +  AfJhi^J(r)] 

n-l 

(2.41) 

Trr  \  =  ^cos{\i„z)Wjh^^Xr)  +  (r)  +  Al^JhgJ(r)] 

rt— 1 

(2.42) 

=  ^cosi[l„z)Wjhlll(r)  +  A^Jhl\l(r)  +  A'i^^hlUir)] 

n—l 

(2.43) 

Tzz\  =  ^cos(\i,z)[Ail,^hlll(r)  +  A^Jhlll(r)  +  A'iJhiZ(r)^ 

n—l 

(2.44) 

Trz  \  =  £  sin(  \i„z)Wnhyl  (r)  +  (r)  +  (r)] 

n—l 

(2.45) 

+  A^y^ir)  +  A^^hlfnir)] 

n—l 

(2.46) 

Dz\  =  ^cos(\l^z)Wjh'^^l(r)+A^2jh'2i(r)  + 

(2.47) 

«=1 


where  ,  and  are  sets  of  constants  and  are  to  be  determined  from  the  boundary 

conditions.  The  |i„ 's  are  the  eigenvalues  and  determined  by  the  boundary  conditions  at  z=±l. 


(n  =  1,3,5,...).  (2.48) 


The  functions  h^l\r)  through  are  given  in  the  Appendix. 


For  the  isotropic  and  nonpiezoelectric  interlayer  (i=2)  and  matrix  (/=3)  the  governing 
equations,  Eqs.  (2.20)  through  (2.22),  reduce  to 


V72,,(0  ^  ,  ^11  ^^12  " 

fdn‘'>  ,  ,  dw^'-'  ^ 

'  ‘'ll 

JO  ,  (0  Pi 

ii  ^r\ 
f  du‘^ 

dr  r 

dw^'^  ^ 

_1_  _1_ 

dz  J 

=  0 

'  ’’  1  /■\  /  ■\ 

c  ^  —  c  ^  3r 

Cfi  Cj2 

dr  r  dz  j 

=  0 


(2.49) 

(2.50) 


u  V72  1  3 

where  V  =  — ^  +  -  —  +  — 
or  rar  az 
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The  complementary  solution  to  Eqs.  (2.49)  and  (2.50)  is  obtained  by  using  Love’s  stress 
function  provided  the  stress  function  F^'\r,z)  satisfies  the  equation, 

V"VV('')  =  0  (/  =  2,3)  (2.51) 

Thus,  the  complementary  solution  for  the  interlayer  (/=2)  and  the  matrix  (/=3)  are  given  by 
uV  =  l^cos(|l„z)[A';;^/z[;’(r)  +  (^)]  (2.52) 

w— 1 

wf  =  £)sin(  |i„ z)Wjh!i'Xr)  +  {r)  +  A^l gl‘J(r)  +  (r)]  (2.53) 

n—\ 

Trr\  =  f^cos(|l„z)[Ai'‘^/z5l’(r)  +  +  A^lgf^ir)  +  Af„gtl{r)\  (2.54) 

=  £cos(|l„z)[A'‘^/z7l’(r)  +  ^^y^'l{r)  +  Af^g\'l{r)  +  Al‘^g^^^(r)]  (2.55) 

n-l 

Tzz\  =  ^cos([l^z)[Al!,^h^g„\r)  +  Al%li(r)  +  A^‘^gl‘^(r)  +  A[‘^gloi(r)]  (2.56) 
Tnc  =  £sin(  [i„z)[Al‘Jhun(r)  + Al%i^Z(r)  +  A^lg^^lir)  +  Atlgyl{r)\  (2.57) 

M— 1 

Tnc  =  I^sin(  \i^z)[A\y\Z{r)  + Af^hlZir)  +  A^lg^Zir)  +  A^^^gl‘l(r)]  (2.58) 

n-l 

where  hZ\r)  through  expressions  containing  the  material  constants  and  the 

modified  Bessel  functions  of  the  first  kind,  /o(|J.„r)  and  /i(|l„r),  and  glZ(f  )  through  glZ(f') 

are  expressions  containing  the  material  constants  and  the  modified  Bessel  functions  of  the 
second  kind,  ^o(|l„r)  and  K^(\x^r).  The  functions  hZ\r)  through  h^2„(r)  and  gZJ(f) 
through  gj|„(r)  (/=2,3)  are  given  in  the  Appendix.  The  constants  A[^\  Aj^,  A^^  and  A^'j 
(i=2,  3)  are  to  be  determined  by  boundary  conditions. 
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The  final  solution  in  each  region  is  achieved  by  superposing  the  complementary  solution  and 
an  additional  solution  corresponding  to  a  generalized  plane- strain  state  in  order  to  satisfy  all  the 
boundary  conditions.  The  generahzed  plane -strain  solution  for  the  ceramic  rod  (/=1)  is  given  by 


(1)., 


(1) 


(1) 

w  z 


(1). 


-d) 


(1) 


T„  g  —  Cn 


^(1)  _]1 


(DA 


+  c 


t(D  _„(0 
-'ee  g  ‘'12 


(D  T1 


J 

(D  A 


(D  A 


(Dz7(D 
z  g 


+  C 


V  '  J 
Dl)e(i)  ,  „(l)r(i) 


^“'+1 


(DA 


^  zz  g  ~  ^‘'13  S  ^  ‘'33  S  <^33 


V 

(1)  I7(D 

z  g 


+  .(Dr  (1)  _e(i)F(D 
^‘'13S  ‘^31 s 


rriyl)  _  rTi{L)  _  rTi{l)  _  7^' _  O 

^  rz  g  -  ^^Z  g  ~  ^  rQ  g~  g  -  ■ 


(2.59) 

(2.60) 
(2.61) 

(2.62) 

(2.63) 

(2.64) 


The  generahzed  plane  strain  solution  for  the  interlayer  (i=2)  and  matrix  (/=3)  is  given  by 


( D  t  ( D  "n 

u  =  q  r  +  — 


(D 


(0  yU) 

w  =q  z 


T-’(0  _  (0 

-^rr  g  ~  ^11 


^  (D  _  (1 


(0  A  r 

+  c^'’ 

^  ‘'12 


T-’(0  _  (0 

''ee  g  ~  ^12 


(0  T1 


J 

(D  A 


(0  A 


V 


(0^  (D 


+  c 


(0 


(0  A 


(Or  (0 


+4^C 


'13 

Y^CO  _  'T'iO  _ Y^IO  _  o 

^  rzg~  ^^.g  r&g~'^- 


(2.65) 

(2.66) 

(2.67) 

(2.68) 

(2.69) 

(2.70) 


The  constants  ^  and  ^  (/=1,  2,  3)  are  to  be  determined  by  the  boundary  conditions. 


•(0 


2.3  Boundary  Conditions 
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The  total  stresses  T.j  =  T^j^+T.j^,  displacements  +  and  the  total  electrical 

potential  (j)  =  required  to  satisfy  the  boundary  conditions.  The  boundary  conditions 

for  the  single  embedded  ceramic  rod  with  interlayer  shown  in  Fig.  2.2  are  as  follows: 

(1)  At  r  =  0 ,  the  solution  is  bounded. 


(2)  At  z  =  ±/ ,  the  normal  traction  is  ,  the  shear  stresses  constitute  a  system  in 

equihbiium  and  the  voltage  apphed  is  ±E^l .  These  conditions  are  written  as 


y^(l)  _  _  y^(3)  _  rpO 

^  zz  ~  ^  zz  ~  ^  zz  ~  ^  zz 


(2.71) 


jr^'Vr^O  (2.72) 

4)=±£'°/.  (2.73) 


(3)  At  r  =  ,  the  normal  surface  traction  is  and  the  tangential  traction  is  free,  i.e. 


y^(3) 

1  ff. 


y^(3) 
^  rz 


=  0. 


(2.74) 

(2.75) 


For  both  the  rod- interphase  boundary  and  the  interphase- matrix  boundary,  continuity  of 
displacements  and  traction  (perfect  adhesion)  is  assumed.  Also  the  ceramic  rod- interphase 
boundary  is  treated  as  a  dielectric -air  boundary.  Thus,  the  interface  boundary  conditions  are 
expressed  as 

(4)  At  r  =  q 


(2.76) 


w 


(2.77) 
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_  0^(2) 

rr  ^  rr 

(2.78) 

0^(1)  _  0^(2) 

^  rz  ~  ^  rz 

(2.79) 

D®  =0. 

rr 

(2.80) 

(5)  At  r=  r2 


(2.81) 

(2.82) 

0^(2)  y^(3) 

rr  ^  rr 

(2.83) 

0^(2)  _  t^(3) 

^  rz  ~  ^  rz  * 

(2.84) 

Application  of  the  boundary  conditions  at  r  =  0  yields 


Tl'«=0 


(2.85) 


(2) 


The  remaining  boundary  eonditions  generate  two  systems  of  hnear  equations  which 

(i=2,  3). 


determine  the  unknown  eonstants  ^ ^ \  ^ 


40;  41 

^In  ’  ^2«  ’  ^3n 


andA;:\  A;„ 


2.4  Static  Performance  of  Piezocomposites:  The  Local  Displacement  Field 

In  general,  the  piezoceramie  and  polymer  in  a  1-3  piezocomposite  wiU  have  different 
deformations  when  the  composite  is  subjected  to  an  external  electrical  or  mechanical  stimulus. 
These  deformations  are  dependent  on  the  relative  eomphance  of  the  ceramie  and  the  polymer 
matrix  and  the  interface  between  them.  The  loeal  deformation  also  determines  the  effective 
behavior  of  the  composite.  The  composite  cyhnder  model  developed  is  used  to  simulate  the 
local  deformations  of  a  1-3  piezoeomposite  driven  by  an  external  electric  field. 
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Consider  a  1-3  piezocomposite  subjected  to  an  electric  field  in  the  form  of 

(2.86) 

(2.87) 


P,(a)=o 

(|,(±/)  =  ±£”/=±| 


where  V  is  the  voltage  applied  to  the  ends  of  the  composite  cy Under.  Using  the  solution 
method  outUned  above,  predictions  were  made  of  the  static  displacements  of  a  single-rod 
piezoelectric  composite  assuming  the  PZT-5H  properties  listed  in  Table  2.1  (Morgan  Matroc 
Inc.,  1989).  Geometric  parameters  are  =0.375  mm,  r2  =0.475  mm,  rj  =  3.0  mm  and 

I  =5.0  mm.  The  voltage  appUed  to  the  ends  of  the  composite  cylinder  is  70  V. 

Table  2.1.  Elastic  and  piezoelectric  properties  of  PZT-5H  rods  . 


Property 

PZT-5H 

Property 

PZT-5H 

Cl?  (GPa) 

126 

c?  (C/m? 

-6.55 

c?  (GPa) 

79.5 

c?  (C/m? 

23.3 

c?  (GPa) 

84.1 

c?  (C/m? 

17.0 

c?  (GPa) 

117 

d?  (10"‘VV) 

-274 

C4?  (GPa) 

23.0 

d?  (io"‘Vv) 

593 

(10''F/m) 

F30 

(io-‘Vv) 

45 

e?'’  (lO'^F/m) 

F50 

83?'^  (lO'^'F/m) 

3.01 
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2.4.1  Comparison  with  simple  models 


The  axial  displacement  at  z  =  /  of  the  piezocomposite  cylinder  is  first  calculated  for  the  case 
of  no  interlayer.  Displacement  values  are  compared  with  those  predicted  by  simple  plane -strain 
models  and  are  plotted  as  a  function  of  radial  distance  in  Fig.  2.3.  The  displacement  profile 
obtained  according  to  the  simple  models  is  uniform  because  all  of  the  simple  models  (Haun  and 
Newnham,  1986;  Chan  and  Unsworth,  1989;  Jensen  1991)  assume  uniform  deformation  in 
each  constituent  and  give  the  same  prediction  of  the  effective  composite  piezoelectric  constant, 
^33 ,  which  is  expressed  as 


^33  —  ■ 


VrS 


(3) 


/033 


-|-(1- V.)5. 


(1) 

33 


(2.88) 


where  is  the  volume  fraction  of  piezoelectric  ceramic  rods  and  are  elastic  compliance 
constants  for  each  phase. 


When  a  1-3  piezocomposite  in  these  simple  models  is  subjected  to  uniform  axial  electric 
field  (£) )  only,  the  axial  strain  is  given  by  and  the  axial  displacement  is  w  =  . 


The  present  solution  predicts  a  nonuniform  displacement  profile.  The  maximum 
displacement  occurs  in  the  rod  and  decreases  rapidly  with  radial  distance  in  the  matrix.  As  the 
length  of  the  rod,  21 ,  approaches  infinity,  the  solution  approaches  that  of  a  generahzed  plane - 
strain  model,  such  as  the  concentric  cylinder  model  by  Jensen  (1991). 

2.4.2  Influence  of  the  interlayer  properties 

A  parametric  study  was  carried  out  to  assess  systematically  the  influence  of  the  interlayer 
properties  on  the  static  displacement  of  the  current  1-3  piezocomposite  model. 
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Displacement,  w  (p  m) 


Radial  distance,  r/r 

1 


Fig.  2.3.  Predictions  using  the  current  model  and  simple  plane- strain  model  of  axial 
displacement  at  z  =  I . 
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The  interlayer  was  chosen  to  have  an  elastic  stiffness  either  higher  or  lower  than  the  neat  matrix. 
The  properties  of  the  various  interlayers  used  in  the  calculations  are  given  in  Table  2.2.  The 
symbol,  in  Table  2.2  represents  the  Young’s  modulus  of  the  interlayer  and  v*^^  is  the 
Poisson’s  ratio.  The  Young’s  modulus  and  the  Poisson's  ratio  of  the  surrounding  matrix  are 
kept  as  =  2.101  GPa  and  v^^^=0. 3.  Four  different  cases  were  considered.  In  Case  1 
the  interlayer  and  the  matrix  are  the  same.  In  Case  2  the  interlayer  is  softer  than  the  matrix  while 
in  Case  3  the  interlayer  is  stiffer  than  the  matrix.  In  Case  4  the  interlayer  modulus  is  one 
hundredth  of  the  matrix  modulus.  The  matrix  is  assumed  to  be  Spurr  epoxy  and  the  properties 
were  calculated  using  data  from  the  manufacturer  and  room  temperature  tension  tests. 

Table  2.2.  Elastic  properties  of  interlayer  used  for  parametric  studies. 


Property 

Case  1 

Case  2 

Case  3 

Case  4 

clf  (CPa) 

2.827 

0.565 

14.135 

0.0287 

(GPa) 

1.211 

0.242 

6.055 

0.01211 

(CPa) 

2.101 

0.420 

10.50 

0.021 

V(2) 

0.3 

0.3 

0.3 

0.3 

Axial  displacement  profiles  at  z  =  I  of  the  piezocomposite  cylinder  corresponding  to 
interlayers  with  various  elastic  properties  are  shown  in  Fig.  2.4.  In  the  calculation,  the 
parameters  r^,  r^,  I  and  V  remain  unchanged.  It  is  seen  in  Fig.  2.4  that  the  axial 

displacement  of  the  PZT  rod  is  significantly  influenced  by  introducing  a  comphant  interlayer  with 
a  thickness  of  only  0.1  mm.  The  axial  displacement  of  the  PZT  rod  increases  as  the  modulus  of 
the  interlayer  decreases.  The  maximum  displacement  of  the  rod  can  be  increased  64%  by 
reducing  the  interlayer  modulus  to  one  hundredth  of  the  matrix  modulus. 
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Displacement,  w  (p  m) 


Radial  distance,  r/r 


Fig.  2.4.  Axial  displacement  profiles  at  z=l  for  interlayers  of  varying  moduli  and 
thickness  of  0.1  mm. 
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2.4.3  Influence  of  the  matrix  properties 


The  influence  of  the  matrix  was  examined  by  keeping  the  interlayer  properties  constant  while 
varying  the  matrix  properties.  Three  different  cases  were  considered.  In  Case  1  the  interlayer 
and  the  matrix  are  the  same.  In  Case  5  the  matrix  is  softer  than  the  interlayer  while  in  Case  6 
the  matrix  b  stiffer  than  the  interlayer.  Properties  of  the  interlayer  are  =  2.101  GPa  and 
=  0.3  .  The  parameters  r^,  r^,  r^,  I  and  V  also  remain  unchanged.  Properties  of  the 

matrix  are  tabulated  in  Table  2.3  and  the  results  are  shown  in  Fig.  2.5.  The  out-of-plane 
displacement  of  the  single-rod  piezocomposite  depends  highly  on  the  modulus  of  the  matrix. 
The  trend  of  displacement  variation  shown  in  Fig.  2.5  agrees  with  that  predicted  by  simple 
models,  i.e.,  a  composite  with  lower  matrix  modulus  has  a  larger  displacement. 

Table  2.3.  Elastic  properties  of  matrix  and  interlayer  used  for  parametric  studies. 


Property 

Case  1 

Case  5 

Case  6 

(GPa) 

2.827 

0.565 

14.135 

c,?  (GPa) 

1.211 

0.242 

6.055 

(GPa) 

2.101 

0.420 

10.50 

v(3) 

0.3 

0.3 

0.3 

(GPa) 

2.101 

2.101 

2.101 

V(2) 

0.3 

0.3 

0.3 

2.4.4  Influence  of  the  interlayer  thickness 

The  study  of  interaction  between  the  polymer  matrix,  the  interphase  and  the  piezoelectric 
ceramic  was  further  developed  by  investigating  the  influence  of  interlayer  thickness  on  the 
displacement  profile.  Thickness  variations  were  studied  for  both  a  soft  interlayer  (Case  2  in 
Table  2.2)  and  a  stiff  interlayer  (Case  3  in  Table  2.2). 
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Displacement,  w  (|am) 


Radial  distance,  r/r 

1 


Fig.  2.5.  Axial  displacement  profiles  at  z=l  for  matrices  of  varying  moduli  (the 
interlayer  thickness  =  0.1  mm). 
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The  thickness  of  the  interlayer  is  defined  as  r^  —  r^.  The  maximum  axial  displacement  which 
is  achieved  at  the  center  of  the  rod  is  plotted  versus  the  interlayer  thickness.  Results  are  shown 
in  Fig.  2.6  and  Fig.  2.7.  For  a  soft  interlayer,  the  displacement  increases  as  the  interlayer 
thickness  increases.  For  a  stiff  interlayer,  the  displacement  decreases  as  the  interlayer  thickness 
increases.  The  results  show  that  the  influence  of  the  interlayer  becomes  more  significant  as  the 
thickness  of  the  interlayer  increases. 


2.4.5  Influence  of  the  ceramic  rod  volume  fraction 

The  volume  fraction  of  ceramic  rods  is  an  important  parameter  for  designing  and  evaluating 

piezocomposites.  In  the  current  single-rod  model,  the  volume  fraction  of  rods  is  calculated  by 
V  f  =  {rj  r^)  .  The  effect  of  increasing  rod  volume  fraction  is  shown  in  Fig.  2.8  where  maximum 

axial  displacement  is  plotted  as  a  function  of  v^.  Interlayer  and  matrix  parameters  used  in  the 
calculation  correspond  to  Case  2  in  Table  2.2.  As  the  volume  fraction  of  rods  increases,  the 
maximum  displacement  also  increases. 


2.5  Effective  Properties  of  Piezocomposites 

Consider  the  overall  behavior  of  the  1-3  piezocomposite  shown  in  Fig.  2.1  exhibiting 
transverse  isotropy.  The  constitutive  relations  between  the  average  field  variables  ( (^y) ,  (^i ) ) 

and  ( (Sij^ ,  (D.))  are  described  by  means  of  the  effective  properties  such  that 
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Max  displacement,  w  (p  m) 


0.0155 


0.0150 


0.0145 


0.0140 


0.0135 


0.0130 


Interlayer  thickness  (mm) 


Fig.  2.6.  Variation  of  the  maximum  axial  rod  displacement  with  interlayer  thickness,  for 
a  soft  interlayer  5). 


33 


Max  displacement,  w  (pm) 


Interlayer  thickness  (mm) 


Fig.  2.7.  Variation  of  the  maximum  axial  rod  displacement  with  interlayer  thickness,  for 
a  stiff  interlayer  (7*^'  =  57^^*). 
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Rod  volume  fraction,  v 
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Fig.  2.8.  Variation  of  the  maximum  axial  rod  displacement  with  volume  fraction  of  PZT 
rods. 
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where  (•)  represents  volume  average  and  an  over  bar  denotes  an  effective  property  for  the 
piezocomposite. 

2.5.1  Effective  piezoelectric  constants  and 

In  order  to  determine  one  can  impose  the  following  boundary  conditions  on  the 
piezocomposite 


Pz(±i)  =  K 

(2.91) 

pXh)  =  ^ 

(2.92) 

4)(a)  =  0 

(2.93) 

which  correspond  to  axial  loading  along  the  x,  axis  of  the  composite  and  no  external  electric 
field  apphed.  The  volume  average  of  the  resulting  dielectric  displacement  in  the  composite 
cylinder  element  shown  in  Fig.  2.2  is  calculated  as 


/ 

\ 


1^(1)  \  _/j  -i-  ^  4-/^  'J^W\  ,  I  ryT 


(2.94) 


The  effective  piezoelectric  constant  of  the  1-3  piezocomposite  is  then  given  by 


^33  =  ^  •  (2.95) 
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To  determine  djp  a  radial  stress  is  applied  to  the  boundary  of  the  piezocomposite  such  that 


0 

II 

(2.96) 

PAh)  =  Tl 

(2.97) 

4)(a)=  0 

(2.98) 

The  volume  average  of  the  dielectric  displacement  in  the  composite  cyhnder  element  is  again  of 
the  form 


-(1)\ 


(2.99) 


The  effective  piezoelectric  constant  of  the  1-3  piezocomposite  is  then  given  by 

dn=^.  (2.100) 

2.5.2  Predictions  of  and 

The  piezoelectric  charge  constants  djj  and  d^^  are  calculated  using  Eqs.  (2.95)  and 
(2.100).  The  constants  d^^  and  d^^  are  plotted  against  the  volume  fraction  of  PZT  rods  in  Fig. 
2.9  for  the  case  of  no  interlayer  and  fixed  PZT  rod  aspect  ratio,  111 d .  Here  d  is  diameter  of 
the  piezoceramic  rod  and  21  is  length  of  the  rod.  Like  the  maximum  static  displacement,  both 
^33  and  -^3;  increase  with  increasing  rod  volume  fraction.  Both  ^33  and  -^3^  are  maximum 

for  the  pure  ceramic  when  the  rod  volume  fraction  is  unity.  However,  some  properties  of  the 
composite,  like  the  effective  hydrostatic  piezoelectric  constant  d^  =  d^^^  +  ld^^  discussed  in 
Chapter  4,  can  be  improved  for  the  intermediate  rod  volume  fraction  region.  The  values  of  ^33 
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and  -^31  predicted  by  the  simple  generalized  plane- strain  model  are  also  plotted  in  Fig.  2.9. 
These  values  are  higher  than  that  obtained  from  the  current  model. 

It  is  shown  in  Fig.  2.10  and  Fig.  2.11  the  influence  of  the  PZT  rod  aspect  ratio  on  the 
effective  piezoelectric  constants  ^33  and  for  the  case  of  no  interlayer.  The  values  of  ^33  and 
-^31  increase  with  increasing  aspect  ratio,  but  there  is  a  saturation  of  the  effect  for  21  Id  >  500 . 
These  saturated  values  of  ^33  and  ^3;  approach  the  values  predicted  by  the  simple  generahzed 

plane- strain  model.  The  influence  of  the  interlayer  modulus  on  the  effective  piezoelectric 
constants  d^^  and  d^^  is  shown  in  Fig.  2.12  and  Fig.  2.13  for  an  interlayer  that  is  0.1  mm  thick. 

The  effect  of  the  interlayer  modulus  depends  on  rod  volume  fraction.  For  an  interlayer  with 
modulus  lower  than  that  of  the  matrix,  d^^  and  -d^^  decrease  for  low  rod  volume  fraction  and 
increase  for  high  rod  volume  fraction  (v^  >  0.5  ).  For  an  interlayer  with  modulus  shghtly  higher 
than  the  matrix  the  d^^  and  -d^^  decrease  for  0  <v^  <0.6  .  For  a  very  comphant  interlayer 
both  ^33  and  -d^^  decrease,  but  the  hydrostatic  piezoelectric  constant  may  increase,  as  will 
be  discussed  in  detail  in  Chapter  4. 

2.6  Discussion 

Difficulties  exist  in  manufacturing  1-3  piezocomposites  with  large  21  Id  ratio.  In  practice, 
values  of  2lld  are  typically  less  than  25.  Thus,  1-3  piezocomposites  should  be  modeled  as 
ahgned  short- fiber  composites  and  not  as  continuous -fiber  composites.  In  this  chapter,  a  three- 
dimensional  analytical  model  was  developed  for  studying  the  local  fields  and  effective  behavior 
in  1-3  piezocomposites. 

The  results  obtained  using  the  current  model  demonstrate  that  the  local  out-of-plane 
displacements  of  the  composite  are  highly  dependent  on  the  interaction  between  the  polymer 
matrix,  a  thin  interlayer,  and  the  piezoelectric  ceramic.  In  order  to  enhance  the  out-of-plane 
displacement,  a  passive  matrix  with  small  Young's  modulus  should  be  chosen.  This  result  is 
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similar  to  that  obtained  from  simple  models.  However,  the  overall  stiffness  of  the  composite  is 
dependent  on  the  matrix  modulus,  especially  for  lower  volume  fractions.  By  introducing  a  thin 
comphant  interlayer,  the  out-of-plane  displacement  of  the  rod  is  greatly  increased  while  the 
requirement  for  overall  stiffness  of  the  composite  can  stiU  be  satisfied.  The  interlayer  essentially 
serves  to  couple  or  decouple  the  rods  from  the  matrix,  depending  on  the  interlayer  stiffness. 
Interlayer  properties  were  also  found  to  influence  significantly  the  effective  charge  constants 
and  djj  for  the  piezocomposite.  Thus,  the  interlayer  can  be  tailored  to  optimize  the  properties 
and  increase  the  sensitivity  of  1-3  piezocomposites. 
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Piezoelectric  constants  m/V) 


Fig.  2.9.  Variation  of  effective  piezoelectric  constants  with  PZT  rod  volume  fraction  for 
the  case  of  no  interlayer. 
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Fig.  2.10.  Variation  of  with  PZT  rod  aspect  ratio  21  Id . 
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Fig.  2.11.  Variation  of  djj  with  PZT  rod  aspect  ratio  2lld . 
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Fig.  2.12.  Variation  of  with  interlayer  modulus  for  an  interlayer  of  0.1  mm 
thickness. 
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Fig.  2.13.  Variation  of  with  interlayer  modulus  for  an  interlayer  of  0.1  mm 
thickness. 
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3.  IN-SITU  DISPLACEMENT  MEASUREMENTS 
OE  PIEZOCOMPOSITES 


3.1  Motivation 

In  Chapter  2,  a  composite  cylinder  model  was  used  to  predict  the  local  deformations  of  a 
1-3  piezocomposite  driven  by  an  external  electric  field.  The  local  deformations  were  highly 
sensitive  to  the  piezoelectric  properties  of  the  rod  and  the  elastic  properties  of  the  polymer.  The 
properties  of  the  interface  region  between  the  rods  and  the  matrix  were  also  shown  to  have  a 
significant  effect  on  local  behavior.  The  in- situ  rod  displacements  induced  by  apphcation  of 
electric  field  (converse  piezoelectric  effect)  are  indicative  of  piezocomposite  performance.  For 
example,  the  average  out- of- plane  surface  displacement  can  be  correlated  with  the  effective 
piezoelectric  constant  for  the  composite. 

In  this  chapter,  the  in- situ  static  displacements  of  a  1-3  piezocomposite  are  investigated 
experimentally  using  laser  interferometry.  Displacement  profiles  are  measured  for  several 
different  interface  conditions  and  are  correlated  with  the  rod/polymer  adhesion.  Experimental 
results  are  also  compared  with  predictions  made  by  using  the  composite  cyhnder  solutions 
derived  in  Chapter  2. 

3.2  Micro -interferometric  Method 

Quantitative  analysis  of  the  local  static  deformations  in  1-3  piezoceramic -polymer 
composites  is  a  challenging  problem.  Since  the  deformations  are  small  and  nonuniform,  the 
experimental  technique  must  have  an  out- of- plane  sensitivity  of  approximately  5  nm  and  a 
spatial  resolution  commensurate  with  the  dimensions  of  the  composite  unit  cell.  Previous 
investigations  by  Sottos,  Scott  and  McCullough  (1991)  and  Ryan,  Scott  and  Sottos  (1990) 
have  demonstrated  the  utihty  of  heterodyne  micro- interferometric  methods  for  nondestmctively 
measuring  displacements  in  the  interfacial  regions  of  polymer  matrix  composites.  Because  a 
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long  distance  microscope  with  a  high  magnification  is  used,  the  interferometer  probe  has  the 
capabihty  to  resolve  small  displacements  in  a  very  small  region.  The  current  investigation  adapts 
this  micro- interferometric  technique  for  measuring  in- situ  surface  displacements  in  1-3 
piezocomposites. 

3.2.1  Apparatus 

A  scanning  heterodyne  micro- interferometer  similar  to  that  described  in  Sottos,  Scott  and 
McCullough  (1991)  is  utihzed  to  measure  the  surface  displacements  of  1-3  piezocomposite 
samples.  A  schematic  diagram  of  the  scanning  heterodyne  micro- interferometer  is  shown  in  Fig. 
3.1.  A  single,  hnearly  polarized  hght  beam  of  wavelength  514.5  nm  from  an  argon  laser  (Lexel 
Laser  model  3500)  is  incident  upon  a  40  MHz  acousto- optic  modulator  (AOM)  producing  two 
beams  which  are  sent  along  different  arms  of  the  interferometer.  The  first  beam,  which  follows 
path  A  in  Fig.  3.1,  is  shifted  in  frequency  by  40  MHz  and  is  used  as  a  reference  beam.  The 
second  beam  with  the  same  frequency  as  the  beam  incident  upon  the  AOM  is  directed  along 
path  B  in  Fig.  3.1  and  is  used  to  illuminate  the  sample  surface.  Hence,  the  sample  serves  as  a 
mirror  in  this  arm  of  the  interferometer.  The  sample  is  mounted  on  a  translation  and  rotation 
stage,  such  that  only  the  outer  edges  of  the  back  face  of  the  sample  are  actually  touching  the 
fixture.  This  arrangement  allows  the  sample  to  displace  freely.  The  sample  is  manually  adjusted 
for  angular  tilt  when  ahgning  the  interferometer. 

The  combination  of  the  polarizing  beam  sphtter  and  the  quarter- wave  plate  allows  the 
sample  to  be  illuminated  at  normal  incidence  without  excessive  hght  loss.  The  polarizing  cube 
reflects  vertically  polarized  light,  and  transmits  horizontahy  polarized  hght.  The  vertical  polarized 
hght  incident  upon  the  cube  is  reflected  through  the  quarter- wave  plate,  which  circularly 
polarizes  it.  After  reflection  from  the  sample,  the  beam  passes  back 
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Fig.  3.1.  Schematic  of  the  scanning  micro- interferometer. 
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through  the  quarter- wave  plate,  producing  a  horizontally  polarized  beam  that  passes  through  the 
polarizing  cube  toward  the  nonpolarizing  beam  sphtter.  The  half- wave  plate  in  the  reference 
arm  rotates  the  polarization  of  the  modulated  beam  to  match  that  of  the  beam  from  the  sample 
arm.  The  two  beams  are  recombined  at  the  nonpolarizing  beam  splitter. 

The  combined  beams  then  arrive  at  the  primary  lens  of  the  long-distance  microscope 
(Infinity  model  K2)  placed  before  the  eyepiece  lens  as  shown  in  Fig.  3.1.  An  image 
(interferogram)  of  the  sample  surface  is  magnified  onto  the  image  plane  with  this  long-distance 
microscope,  where  a  scanning  photodiode  detector  scans  the  magnified  image.  The  large 
numerical  aperture  of  the  long  distance  microscope  increases  light  gathering  power  and  permits 
high- resolution  imaging.  Through  the  use  of  different  lenses,  magnifications  from  lx  to  lOOx  are 
attainable.  Depending  on  the  magnification  of  the  image,  the  in- plane  resolution  is  limited  either 
by  the  optical  system  or  by  the  diameter  of  the  detector  window. 

A  dynamic  image  which  travels  at  a  40  MHz  beat  signal  across  the  image  plane  is  formed 
because  of  the  difference  in  frequency  between  the  sample  and  reference  beams.  One  comer  of 
the  image  plane  is  reflected  by  a  mirror  onto  a  second  photodetector.  The  signal  at  this 
stationary  photodetector  is  used  as  a  phase  reference  signal.  The  phase  of  the  image  signal  is 
continuously  compared  with  the  phase  of  the  leference  signal.  As  a  result,  such  dynamic 
disturbances  as  air  turbulence,  mechanical  vibration,  and  thermal  drift,  that  occur  within  the 
interferometer  but  away  from  the  sample,  do  not  affect  the  measurements. 

The  signal  processing  unit  of  the  interferometer  is  shown  in  Fig.  3.2.  The  scanning 
photodiode  detector  is  moved  across  the  magnified  sample  image  using  two  actuators 
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Fig.  3.2.  Signal  processing  unit  of  the  scanning  interferometer. 
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which  are  perpendicular  to  each  other.  The  actuators  (Newport  model  850B)  are  controlled  by 
a  microcomputer- controlled  motion  controller  (Newport  model  PMC200P).  A  high-frequency 
lock-in  amplifier  (EG&G  PARC  model  5202)  is  used  to  measure  the  phase  of  the  image  signal. 
The  output  signal  from  the  lock-in  amphfier  is  digitized  by  a  digital  oscilloscope  (Tektronix 
model  TDS  420).  Both  the  oscilloscope  and  the  motion  controller  are  interfaced  with  a  central 
control  unit,  a  Power  Macintosh  (Apple  model  8100/80)  via  the  IEEE-488.2  interface.  The 
sweep  rate  of  the  oscilloscope  is  synchronized  to  the  scanning  of  the  photodetector  by  a  Virtual 
Instmment  (VI)  written  in  EabVIEW  (National  Instmments)  graphic  language. 

3.2.2  Principle  of  operation 

Eor  the  calculation  of  the  change  of  optical  path  length  produced  by  out-of-plane  motions, 
incidence  wave  is  assumed  to  be  plane,  which  is  a  reasonable  assumption  in  the  focal  zone  a 
Gaussian  beam.  With  reference  to  Pig.  3.3,  the  formula  which  relates  the  change  of  optical  path 
length  5  A  to  surface  displacement  u  can  be  derived  as  (Monchalin  et  al.,  1989) 


5A=(k'-k”)*u  (3.1) 

Here  k'  is  the  unit  vector  in  the  incident  direction  while  k°  is  that  of  the  observation  direction. 
Eor  detecting  out-of-plane  displacement  w  ,  as  in  this  study,  the  angle  between  k'  and  k"  is 
180°.  Thus  Eq.  (3.1)  is  simphfied  to 

5A=2w  (3.2) 
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Let  denote  the  sample  beam  and  A^  be  the  reference  beam.  On  the  image  plane  (xy 
plane)  A^  and  A^  are  expressed  as 

4  =  exp  \i\2nfj  +  \\f^  (x,  y)]] 

A  =  y)  exp  |i'[27t  f^t  +  v|/  ^(x,  y)]] 

where  a^(x,y)  and  a^(x,y)  are  the  amphtudes  of  A^  and  A^  respectively. 

The  intensity  of  the  superposed  sample  wave  and  reference  wave,  which  is  detected  by  the 
photodetector,  is  given  as 

I  =  al  +al  +2a^a^cos[2n(f^  -  /o)^ +  (¥r  “¥  j]-  (3-5) 

Here  the  frequency  difference  is  set  by  the  AOM  at  40  MHz.  The  phase  of  the 

scanning  detector  output  current,  \|/  =  \|/^  -  \|/^ ,  is  measured  with  reference  to  the  stationary 

detector.  Any  out-of- plane  displacement  of  the  sample  surface  produces  a  phase  shift  in  the  40 
MHz  beat  signal  at  the  corresponding  image  plane.  The  relation  between  this  phase  shift  At]/ 
and  the  change  of  optical  path  length  5  A  due  to  out- of- plane  displacement  is 

(3.6) 

where  X  is  the  wavelength  of  the  laser  hght.  The  displacement  w  can  then  be  calculated  from 
this  phase  shift.  Obtaining  a  A-/100  =  5  nm  sensitivity  would  require  a  phase  measurement 
accuracy  of  approximately  7°,  which  is  guaranteed  by  the  lock-in  amphfier  used. 

3.3  Sample  Preparation  and  Characterization 

3.3.1  Sample  fabrication 


(3.3) 

(3.4) 
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The  interferometer  is  used  to  measure  the  static  deformations  resulting  from  activating  a 
single  PZT  rod  embedded  in  a  polymer  matrix.  A  schematic  of  the  specimen  is  shown  in  Fig. 
3.4.  The  piezoceramic  (PZT)  rods  produced  by  Ceramic  Processing  Systems  were  1  mm  in 
diameter  and  poled  in  an  oil  bath  under  an  electric  field  of  20  kV/cm.  The  axial  piezoelectric 
constant  was  measured  as  =  420  x  10“^^  m/V  .  For  the  matrix,  a  Shell  EPON-828  epoxy 

resin  with  a  Shell  DETA  curing  agent  was  chosen.  The  samples  were  fabricated  by  placing  a 
PZT  rod  into  the  center  of  a  mold  which  was  then  filled  with  epoxy.  The  samples  were  cured  at 
room  temperature  for  7  days  and  post  cured  at  40  “C  for  4  hours  to  raise  the  glass  transition 
temperature  of  the  epoxy. 

The  influence  of  a  distinct  interface  region  between  the  rod  and  the  polymer  was 
investigated  by  fabricating  samples  with  several  different  interface  conditions.  In  the  first  type 
(I),  the  rods  were  embedded  with  no  special  interface  treatment  or  coating  as  described  above. 
In  the  second  type  (11),  a  thin  polymeric  coating  with  a  Young’s  modulus  lower  than  that  of  the 
bulk  matrix  was  applied  to  the  rods  prior  to  being  embedded.  Type  IH  samples  had  an  even 
lower  modulus  (almost  mbber-hke)  coating  applied  to  the  rods.  PinaUy  for  the  fourth  sample 
type  (IV),  the  rods  were  treated  with  a  silane  coupling  agent  obtained  from  Tiber  Materials  Inc. 

Spurr  epoxy  (Polysciences)  was  used  for  the  compliant  coatings  (types  II  and  HI).  Spurr 
epoxy  is  a  four- component  system.  The  modulus  was  varied  by  changing  the  concentration  of 
the  fiexibilizer--diglycidyl  ether  of  polypropylenglycol  (DER)  in  the  stoichiometric  mix.  Table 
3.1  lists  the  various  compositions  of  Spurr  used  in  this  study  and  the  designation  of  each 
composition.  Coatings  were  applied  by  dipping  the  rods  in  the  epoxy  and  curing  at  70  “C  for  8 
hours.  To  increase  the  thickness,  a  second  layer  of 
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Fig.  3.4.  Schematic  of  the  sample  for  interferometric  measurements. 
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coating  was  applied  to  the  first  layer.  Properties  of  the  coatings  as  measured  by  uniaxial  tension 
tests  are  hsted  in  Table  3.2. 


Table  3.1. 

Spurr  epoxy  compositions. 

Coating 

VCD 

DER 

NSA 

DMAE 

Type  I 

10  g 

6g 

26  g 

0.4  g 

Type  n 

10  g 

9.1  g 

26  g 

0.4  g 

Type  in 

10  g 

19.6  g 

26  g 

0.4  g 

Table  3.2.  Properties  of  the  Spurr  epoxy  coatings. 

Coating 

Young’s  modulus  (GPa) 

Type  I 

2.1 

Type  n 

0.6 

Type  in 

0.01 
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3.3.2  Interfacial  adhesion 


It  is  essential  to  have  some  measure  of  the  interfacial  bond  strength  between  the  coated  PZT 
rods  and  the  polymer  matrix  for  the  evaluation  of  mechanical  response  and  the  development  of 
weU- designed  interfaces.  Displacement  measurements  are  sensitive  to  the  PZT/polymer 
adhesion  as  well  as  the  elastic  properties  of  the  coating.  Although  the  hterature  abounds  with 
studies  of  fiber/mattix  adhesion  in  polymer  composites,  very  httle  work  has  been  carried  out  on 
piezoceramic/polymer  adhesion. 


A  thin- section  push- out  test  was  adopted  to  measure  the  PZT  rod/polymer  interface 
strength  (Li,  Durkin  and  Sottos,  1994).  This  test,  which  has  been  used  extensively  to  measure 
the  in- situ  interfacial  shear  strength  in  fibrous  composites,  is  useful  for  evaluating  the  effects  of 
coatings  on  debonding.  A  schematic  of  the  thin- section  push- out  apparatus  is  shown  in  Fig. 
3.5.  Single  rod  samples  as  described  in  Section  3.3.1  were  cut  to  a  length  of  approximately  3.5 
mm.  Tests  were  conducted  by  applying  an  axial  load  to  the  rod  using  the  punch,  which  was 
attached  to  an  Instron  testing  machine.  A  constant  cross- head  displacement  rate  of  0.05 
mm/sec  was  used  for  the  tests.  A  typical  load  versus  displacement  curve  is  shown  in  Fig.  3.6. 
The  maximum  load  corresponds  to  the  point  at  which  the  rod  debonded  from  the  matrix.  Using 
the  maximum  load,  ,  and  the  sample  dimensions,  one  can  calculate  the  average  interfacial 
shear  strength,  Xg ,  expressed  as 


47lr/ 


(3.7) 


Although  this  simplified  equation  only  calculates  an  average  interfacial  strength,  it  is  useful  for 
comparing  adhesion  changes  due  to  the  different  interlayers  or  surface  treatments  used  in  this 
study.  The  effects  of  different  coating  types  on  interfacial  shear  strength  are  summarized  in 
Table  3.3.  The  results  show  that  the  addition  of  a  interlayer  with  lower  elastic  modulus  has 
almost  no  effect  on  the  average  interfacial  shear  strength 
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Fig.  3.5.  Schematic  of  the  thin- section  push-out  test. 
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Fig.  3.6.  A  typical  load- displacement  curve  for  the  push-out  test. 
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at  room  temperature.  The  results  also  show  that  the  interfacial  shear  strength  is  significantly 
increased  by  the  silane  coupling  agent. 


Table  3.3  Average  interfacial  shear  strength. 


Coating 

Interfacial  Shear  Strength 

Type  I 

17.3  MPa 

Type  II  &  Type  IH 

18.7  MPa 

Type  rV 

22.4  MPa 

3.3.3  Surface  preparation  for  interferometric  measurement 

For  the  interferometric  experiments,  single  rod  samples  were  cut  to  a  length  of  5  mm  and 
the  front  face  finely  pohshed.  Polishing  of  the  face  on  which  the  laser  probe  measurements  are 
made  was  critical  to  insure  a  fiat  surface  free  of  scratches.  After  pohshing,  a  thin  gold  coating 
was  evaporated  on  the  front  surface  of  the  sample  to  provide  a  highly  reflective  surface  for  the 
interferometric  measurement  as  well  as  an  electrode  for  applying  the  electric  field.  Silver  paint 
was  applied  to  the  back  surface  of  the  sample  to  provide  a  second  electrode. 

3.4  Experimental  Procedure  and  Results 

The  piezoceramic  rod  was  actuated  by  applying  a  static  voltage  across  the  piezocomposite 
sample.  Two  types  of  displacement  measurements  are  reported  below.  In  the  first  type  of 
measurement,  maximum  displacement  is  recorded  as  a  function  of  time  after  actuation.  In  the 
second  type  of  measurement,  displacement  profiles  of  the  single-rod  composite  are  obtained. 
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3.4.1  Maximum  surface  displacement 

For  the  maximum  displacement  measurement,  the  scanning  detector  is  fixed  in  the  center  of 
the  PZT  rod.  An  electric  field  is  apphed  across  the  sample,  which  triggers  the  acquisition  of  the 
phase  data  as  the  sample  responds.  The  phase  change  is  then  converted  to  displacement  for 
that  particular  point.  Measured  displacements  are  plotted  as  a  function  of  time  for  aU  four  types 
of  sample  in  Fig.  3.7.  The  apphed  voltage  is  also  shown  in  this  figure.  The  composite  cylinder 
model  discussed  in  Chapter  2  predicts  that  the  maximum  displacement  wiU  occur  at  the  center 
of  the  ceramic  rod.  It  can  be  seen  in  Fig.  3.7  that  the  displacement  at  the  center  of  the  PZT  rod 
jumps  rapidly  in  about  600  ms  from  zero  to  its  stable  value  when  the  dc  voltage  is  turned  on. 
The  displacements  of  the  samples  with  a  low- modulus  interphase  (types  n  and  HI)  were  higher 
than  that  of  the  untreated  samples  (type  I).  Because  the  adhesion  between  the  Spurr  epoxy 
coated  rods  is  very  similar  to  that  of  the  untreated  rods  (see  Table  3.3),  the  increase  in 
displacement  is  due  solely  to  the  reduction  in  stiffness  of  the  interlayer  surrounding  the  rod.  The 
samples  treated  with  a  silane  coupling  agent  (type  IV)  had  a  lower  displacement  than  the 
untreated  samples  (type  I).  This  decrease  in  displacement  may  be  attributed  to  the  increase  in 
adhesion  between  the  silane  treated  rod  and  the  surrounding  matrix  or  to  the  formation  of  a 
stiffer  interlayer  near  the  rod  due  to  a  reaction  between  the  silane  coupling  agent  and  the  matrix 
during  sample  cure.  The  extent  of  reaction  of  the  silane  with  the  matrix  is  not  well  known  and 
therefore  makes  it  very  difficult  to  assess  the  properties  near  the  rod  surface. 

3.4.2  Surface  displacement  profile 

Surface  displacement  profiles  are  obtained  by  scanning  the  detector  in  a  radial  path  from  the 
rod  center.  To  measure  the  net  displacements  of  the  sample,  one  makes  an  initial  scan  of  the 
surface  before  activation.  The  sample  is  then  activated  and  a  second  scan  is  mn.  During  each 
scan,  phase  data  are  continuously  recorded  and  stored.  Figure  3.8  is  a  typical  plot  of  the  raw 
data  from  the  surface  profiles  before  and  after  activation.  Displacements  are  calculated  by 
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subtracting  the  phase  data  of  the  initial  scan  lines  from  the  final  ones.  AH  out-of-plane 
displacements  are  measured  relative  to  the  same  reference  point,  which  is  at  a  distance  at  least  3 
times  the  rod  radius  away  from  the  rod  center. 

Figure  3.9  shows  the  measured  displacement  profiles  along  the  radial  direction  of  the 
samples  for  each  interphase  type.  A  nonuniform  experimental  displacement  profile  is  observed 
with  maximum  displacement  in  ihe  PZT  rod  region.  The  surface  displacement  decreases  with 
radial  distance  in  the  matrix.  The  axial  displacement  of  the  PZT  rod  is  significantly  influenced  by 
a  thin  comphant  coating.  The  rod  displacement  increases  as  the  modulus  of  the  coating 
decreases.  This  observation  is  consistent  with  the  results  predicted  theoretically.  Although  a 
thin  comphant  coating  decouples  the  rod  from  the  matrix  and  increases  the  rod  displacement,  it 
also  increases  the  nonuniformity  of  the  displacement  profile  across  the  composite.  On  the  other 
hand,  samples  with  silane  coupling  agent  have  smaller  and  more  uniform  surface  displacements. 

3.5  Discussion 

A  one  to  one  comparison  of  theory  and  experiment  is  difficult  because  the  mechanical 
properties  of  the  interphase  cured  with  silane  couphng  agent  and  the  thicknesses  of  the  coatings 
are  unknown.  However,  trends  in  the  variation  of  surface  displacement  profile  due  to  changes 
in  interphase  properties  can  be  assessed.  Predictions  of  surface  displacement  profile  were 
made  by  using  the  theoretical  analysis  derived  in  Chapter  2.  Elastic  properties  for  the  PZT  rod 
in  Table  2.1  are  used.  The  measured  value  of  =  420  x  10“*^  m/V  is  used  and  the  value  of 
Jjj  is  estimated  to  be  =  -198  x  10“*^  m/V.  Material  properties  for  EPON  828  epoxy  are 
fisted  in  Table  3.4.  Geometric  parameters  are  =0.5  mm  ,  rj  =0.6  mm  ,  rj  =  3.0  mm  , 
and  I  =2.5  mm  .  The  displacement  profiles  for  samples  with  different  interphase  modulus  are 
shown  in  Pig.  3.10  for  an  applied  dc  voltage  of  300V.  These  displacements  are  relative  to, 
r  =  3rj ,  the  point  where  the  reference  detector  is  located  in  the  experiments.  The  variations  of 
surface  displacement  profile  due  to  changes  in  interphase  properties  as  shown  in  Pig.  3.10  are  in 
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good  agreement  with  the  experimentally  observed  trends  for  the  untreated  samples  (type  I)  and 
samples  with  a  low  modulus  interlayer  (types  n  and  HI).  The  measured  displacement  profile  for 
the  silane  coupling  agent  treated  samples  (type  IV)  is  lower  than  the  predicted  displacement 
profile  for  samples  with  a  high  modulus  interlayer  as  shown  in  Fig.  3.10.  This  decrease  in  the 
measured  displacement  may  be  attributed  to  the  increase  in  adhesion  between  the  silane  treated 
rod  and  the  surrounding  matrix. 

Table  3.4.  Materials  properties  for  EPON  828  epoxy. 


(GPa) 

v(3) 

2.5 

0.3 

Micro- interferometric  measurements  as  well  as  the  theoretical  predictions  demonstrate  that 
the  out- of- plane  displacements  of  the  composite  are  highly  dependent  on  the  interaction 
between  the  polymer  matrix,  a  thin  interlayer  and  the  piezoelectric  ceramic.  By  introducing  a 
thin  comphant  interlayer,  the  out- of- plane  displacement  of  the  rod  is  greatly  increased  while  the 
requirement  for  overall  stiffness  of  the  composite  can  stiU  be  satisfied.  The  interlayer  essentially 
serves  to  couple  or  decouple  the  rods  from  the  matrix,  depending  on  the  interlayer  .stiffness.  On 
the  other  hand  a  silane  coupling  agent  can  improve  the  interfacial  adhesion  and  reduce  the 
nonuniformity  of  the  displacement  profile  in  whole  composite  region.  The  measured 
displacement  profile  can  also  be  used  to  estimate  the  effective  piezoelectric  constants  of  the 
composite.  For  thin- section  1-3  piezocomposites,  can  be  calculated  from  the  out-of- plane 
displacement  as 
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V 


(3.8) 


where  V  is  the  voltage  applied  and  the  integral  is  over  the  area,  Q. ,  measured  on  the  surface  of 
the  composite  sample. 
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Displacement  (pm) 


Time  (ms) 


Fig.  3.7.  Displacements  measured  at  the  center  of  the  PZT  rod. 
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Voltage  (v) 


Surface  profile  (pm) 
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Fig.  3.8.  Measured  composite  sample  profiles. 
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Fig.  3.9.  Displacement  profiles  measured  on  the  composite  sample  surface. 


66 


Surface  displacement  (pm) 


0.1000 


0.0800 


0.0600 


0.0400 


0.0200 


0.0000 


Radial  distance,  r/r 

1 


Fig.  3.10.  Displacement  profiles  predicted  by  theoretical  analysis  (relative  to  r  =  3rj  ). 
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4.  ELECTROMECHANICAL  COUPLING  IN  PIEZOCOMPOSITES 


4.1  Introduction 

The  key  principle  in  designing  piezoelectric  ceramic -polymer  composites  is  to  maximize  the 
electromechanical  coupling.  Maximum  coupling  is  achieved  by  building  a  composite  stmcture 
that  transfers  the  applied  external  stresses  to  the  composite’s  piezoelectric  component  in  a 
pattern  that  most  nearly  approximates  the  piezoceramic’s  maximally  coupled  stress  patterns 
(Smith,  1992).  In  this  chapter,  this  principle  is  apphed  to  the  hydrostatic  behavior  of  1-3 
piezocomposites. 

Piezoceramics,  which  are  used  commonly  as  ultrasonic  transducer  materials,  suffer  from 
several  disadvantages  when  used  for  hydrophone  apphcation.  The  piezoelectric  response  of  a 
PZT  ceramic  under  hydrostatic  pressure  is  shown  schematically  in  Fig.  4.1.  The  stress  induced 
electric  displacement  Dj  can  be  resolved  into  two  parts.  The  first  part  corresponds  to  the  axial 

stress  while  the  second  part  corresponds  to  the  lateral  stress.  Due  to  the  opposite  sign  of  the 
piezoelectric  charge  coefficients  ^33  and  these  two  parts  are  opposite  in  sign  and  the  net 
piezoelectric  effect  is  small.  The  hydrostatic  piezoelectric  coefficient  of  PZT  is  defined  as 
df^  =  2^3  j  -I-  ^33 .  The  coefficient  d^  of  PZT  is  low  even  though  the  magnitudes  of  both  ^33  and 
^3;  are  large.  The  hydrostatic  voltage  coefficient  is  also  small  because  of  the  high 

833 

dielectric  constant  83^3 . 

Smith  (1992)  has  shown  that  the  combination  of  a  lateral  stress  with  sign  opposite  to  that  of 
the  axial  stress  provides  the  maximal  couphng  to  in  a  PZT  ceramic  as  shown  in  Fig.  4.2. 
Thus,  an  effective  way  to  increase  the  magnitude  of  and  gf^  is  the  use  of  a  composite 
stmcture  design.  Among  various  types  of  composite,  piezoelectric  composites  with  1-3 
connectivity  have  attracted  a  great  deal  of  attention  and  have  been  used  widely  for  underwater 
hydrophone  apphcations  (Haun  and  Newnham,  1986). 
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Fig.  4.1.  Schematic  of  the  hydrostatic  response  of  a  piezoceramic. 
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Fig.  4.2.  Schematic  of  the  maximum  direct  stress  pattern  in  a  PZT. 
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The  hydrostatic  response  of  a  1-3  piezocomposite  is  shown  schematically  in  Fig.  4.3.  The 
volume  average  of  the  stress- induced  electric  displacement  (Dj)  is  dependent  on  the  volume 

average  axial  stress  and  volume  average  lateral  stress  in  the  piezoceramic  rods.  Because  of  the 
stress  transfer  from  the  piezoelectrically  passive  polymer  phase  to  the  piezoelectrically  active 
ceramic  phase,  the  hydrostatic  response  of  this  piezocomposite  can  be  dramatically  improved 
over  the  single-phase  piezoceramic.  In  designing  1-3  piezocomposites,  the  goal  of  maximizing 
the  electromechanical  coupling  can  be  achieved  either  by  increasing  the  axial  stress  or  by 
effectively  reducing  the  lateral  stress  in  the  piezoceramic  rods. 

Maximizing  electromechanical  coupling  in  piezocomposites  is  an  inverse  problem.  In  the 
following  two  sections,  the  electromechanical  couphng  in  1-3  piezocomposites  is  quantified  in 
terms  of  the  effective  hydrostatic  piezoelectric  constants  and  .  Since  there  is  no  universal 

method  for  solving  inverse  problems,  a  parametric  study  is  performed  to  examine  systematically 
the  influence  of  matrix  .stiffness,  rod  aspect  ratio,  interlayer  .stiffness,  and  rod  volume  fraction  on 
the  hydrostatic  performance  of  a  1-3  piezocomposite  in  terms  of  the  effective  hydrostatic 
piezoelectric  coefficients. 

4.2  Effective  Hydrostatic  Performance  of  1-3  Piezocomposites 

To  study  electromechanical  couphng  in  1-3  piezocomposites  for  hydrostatic  apphcations, 
the  foUowing  boundary  conditions  are  imposed  on  the  piezocomposite  cyhnder  developed  in 
Chapter  2: 


p^(±/)  =-p 

pXh)  =  -p  ■ 

(^(a)  =  0  , 

where  p  represents  constant  pressure. 


(4.1) 
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Fig.  4.3.  Schematic  of  the  hydrostatic  response  of  a  1-3  piezocomposite. 
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The  volume  average  of  the  hydrostatic  pressure- induced  electric  displacement  in  the 


composite  cylinder  is  given  by 


I  J  7’(l)\  I  /  j  7<1)\  I  /  j  t’CON  I  I  pCON 


(4.2) 


where  (•)  represents  volume  average.  For  a  homogeneous  PZT  rod  with  a  hydrostatic 
piezoelectric  constant  ,  the  volume  average  of  the  pure  hydrostatic  pressure- induced  electric 

displacement  is 


A  =  -^hP  ■ 


(4.3) 


Thus  the  effective  hydrostatic  piezoelectric  constant  of  the  1-3  piezocomposite  is  given  by 


d. 


(4.4) 


The  effective  dielectric  constant  833  of  the  composite  can  also  be  predicted  by  using  the 
composite  cyhnder  model.  Since  the  dielectric  constant  of  the  polymer  is  very  small  when 
compared  with  that  of  the  piezoceramic,  the  effective  dielectric  constant  is  approximated  as 

83^-8[;‘v,  (4.5) 


j'  ( 1) 

where  833  is  the  dielectric  constant  of  the  piezoceramic  and  is  the  volume  fraction  of  the 
piezoceramic  rods.  The  effective  piezoelectric  voltage  constant  is  given  by 


(4.6) 
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The  composite  figure  of  merit  defined  as  is  approximated  as 


(4.7) 

833  V  ^ 

Predictions  were  made  of  the  effective  piezoelectric  constants  and  the  figure  of  merit  in  1-3 
piezoelectric  composites  using  the  Eqs.  (4.4),  (4.6)  and  (4.7).  The  piezoceramic  rods  are 
assumed  to  be  PZT-5H  (Morgan  Matroc  Inc.,  1989)  and  the  properties  are  listed  in  Table  2.1. 
The  polymer  matrix  is  assumed  to  be  Spurr  epoxy  with  the  properties  which  are  given  in  Table 

2.3  (Case  1).  The  diameter  of  the  PZT-5H  rods  is  0.75  mm. 

4.3  Load  Transfer  and  Stress  Amplification 

4.3.1  Predictions  of  effective  dj^  and  for  a  1-3  piezocomposite 

Predictions  of  the  effective  hydrostatic  piezoelectric  constant  for  a  typical  PZT-5H  rod 
aspect  ratio  (21  jd  =  13.3  )  using  the  finite  concentric  cylinder  model  are  plotted  in  Fig.  4.4  for  a 

range  of  rod  volume  fractions  (no  interlayer).  Here  21  is  the  length  of  the  PZT  rods  and  d  is 
the  diameter.  The  effective  exceeds  the  of  the  constituent  ceramic,  approaching  a  value 
three  times  as  high  as  that  of  pure  PZT-5H  (d^  =  45  x  10“'^  m/ V).  The  prediction  of  d^  made 
from  the  generalized  plane-strain  model  is  also  shown  in  Fig.  4.4.  The  value  corresponding  to 
the  generahzed  plane- strain  model  is  larger  than  that  of  the  finite  concentric  cylinder  model. 
Additionally  the  peak  value  of  d^  shifts  to  a  lower  volume  fraction  region  for  the  plane -strain 

prediction. 

Calculated  values  of  and  are  plotted  in  Fig.  4.5  and  Fig.  4.6,  respectively,  as  a 
function  of  the  rod  volume  fraction  for  2/ /d  =13.3.  The  g^  peaks  at  low  ceramic 
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Fig.  4.5.  Variation  of  gf,  with  volume  fraction  of  PZT  rods  for  the  case  of  no  interlayer. 
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Fig.  4.6.  Variation  of  with  volume  fraction  of  PZT  rods  for  the  case  of  no 
interlayer. 
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content,  attaining  a  value  greater  than  that  for  the  pure  ceramic  =  1. 5  x  10“^  Vm/ N).  This 
enhancement  is  due  to  the  dilution  of  the  dielectric  permittivity  and  the  increase  of  d^.  A 
somewhat  broader  peak  in  the  figure  of  merit,  ,  also  appears  at  low  PZT-5H  content  for 
essentially  the  same  reason  as  its  g^  peak.  The  value  attained  is  significantly  higher  than  the 
pure  ceramic  (d^g^  =  67. 5  x  10”^^  m  VN).  Again,  the  plane  strain  model  predicts  significantly 
higher  values  of  and  than  the  composite  cylinder  model. 

The  stress  pattern  in  the  piezoceramic  is  depicted  in  terms  of  the  normalized  average 
stresses  in  the  PZT  rods  which  are  plotted  in  Fig.  4.7  along  with  the  for  the  same  case  as  in 

Fig.  4.4.  There  is  a  large  stress  amphfication  in  the  axial  stress  especially  for  the  low  PZT 
rod  volume  fraction  region.  Although  the  lateral  stress  is  also  shghtly  increased,  competition 
among  the  axial  stress,  the  lateral  stress  and  the  volume  fraction  of  the  PZT  rod  gives  an 
enhanced  effective  hydrostatic  piezoelectric  constant  .  Thus,  axial  stress  amphfication  is 
beneficial  for  the  enhancement  of  . 

4.3.2  Influence  of  PZT  rod  aspect  ratio  on  d^ 

The  aspect  ratio  of  the  PZT  rods,  21  Id ,  has  been  proven  in  experiments  to  be  a  critical 
parameter  in  1-3  piezocomposites.  The  effect  of  the  aspect  ratio  is  not  included  in  the  simple 
generahzed  plane-strain  models  (Smith,  1993;  Jensen,  1991;  Cao,  Zhang  and  Cross,  1992). 
However,  this  effect  can  be  studied  using  the  current  finite  concentric  cylinder  model.  A 
parametric  study  was  carried  out  to  assess  the  influence  of  PZT-5H  rod  aspect  ratio  on  d^ . 

Material  properties  of  PZT-5H  and  the  polymer  matrix  are  the  same  as  in  Table  2.1  and  Table 
2.3  (Case  1).  The  results  are  plotted  in  Fig.  4.8  for  several  different  aspect  ratios.  For  each 
case,  the  exhibits  a  peak  value  which  shifts  to  higher  volume  fractions  with  decreasing  aspect 

ratio.  The  peak  value  increases  with  increasing  aspect  ratio,  but  there  is  a  saturation  of  the 
effect  for  lljd  >  500  .  Thus,  the  enhancement  of  axial  stress  in  the  PZT  rod  from  increasing  the 
rod  aspect  ratio  has  an  upper  limit.  This 
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Fig.  4.7.  Variation  of  d,^  and  the  average  stresses  in  the  PZT  rod  with  rod  volume 
fraction  for  the  case  of  no  interlayer. 
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Fig.  4.8.  Variation  of  d,^  with  volume  fraction  of  PZT  rods  for  different  rod  aspect 
ratios. 
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saturated  value  of  approaches  the  value  predicted  by  the  generahzed  plane  strain  model. 

For  a  composite  with  small  rod  aspect  ratio,  the  simple  generalized  plane  strain  models 
significantly  over  predict  the  .  In  practical  applications,  the  PZT  rod  aspect  ratios  are 

typically  small  {21  jd  <  25)  due  to  poling  limitations,  and  the  finite  concentric  cylinder  model 
would  more  accurately  predict  the  behavior. 

4.3.3  Influence  of  matrix  stiffness  on  d^. 

The  influence  of  matrix  stiffness  on  of  1-3  piezocomposites  was  examined  by  varying  the 
matrix  stiffness  (11  Id  =13.3 )  and  using  the  finite  concentric  cylinder  model.  The  results  are 
plotted  in  Fig.  4.9  where  is  the  Young’s  modulus  of  the  matrix  and  T°  is  the  modulus  of 
plain  Spurr  epoxy  fisted  in  Table  2.3  (Case  1).  The  d^  increases  with  decreasing  matrix 

stiffness  but  saturates  for  a  very  soft  matrix.  The  effect  of  softening  the  matrix  is  similar  to  the 
effect  of  increasing  the  PZT  rod  aspect  ratio.  These  two  effects  are  responsible  for  causing 
axial  stress  amplification,  which  is  significant  at  low  PZT  rod  volume  fractions. 

4.4  Stress  Reduction  and  The  Role  of  Interphase 

The  enhancement  of  electromechanical  coupling  under  hydrostatic  pressure  is  achieved 
either  by  increasing  the  axial  stress  or  by  decreasing  the  lateral  stress  in  the  piezoceramic  rods. 
The  aforementioned  results  of  the  influence  of  the  rod  aspect  ratio  and  the  matrix  .stiffness  on  d^ 

demonstrates  how  axial  stress  amplification  in  the  rods  can  be  achieved  from  these  two  effects. 
The  results  also  show  the  saturation  from  these  two  effects. 

Because  the  actual  stress  transfer  is  through  the  interface  between  the  rods  and  the  matrix, 
the  incorporation  of  any  existing  interphase  region  into  the  micromechanical  analysis  of  1-3 
piezocomposites  may  be  important  for  understanding  composite  behavior. 
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Fig.  4.9.  Variation  of  d,^  with  volume  fraction  of  PZT  rods  for  different  matrix  moduli. 
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Experiments  in  Chapter  3  verified  that  an  interlayer  significantly  influences  the  out- of- plane 
displacement  of  a  1-3  piezocomposite.  In  Chapter  2,  a  thin  interlayer  was  introduced  between 
the  piezoceramic  rods  and  the  matrix  in  the  micromechanical  model.  Interlayer  properties  had  a 
significant  effect  on  the  effective  charge  constants,  ^33  and  ^3^,  as  shown  in  Figs.  2.12  and 

2.13.  It  is  proposed  that  the  presence  of  this  tiin  interlayer  can  influence  the  hydrostatic 
behavior  of  the  composite  and  change  the  sensitivity.  Kim,  Rittenmyer  and  Kahn  (1993) 
showed  that  by  casting  a  soft  epoxy  interlayer  between  PZT  rods  and  a  stiff  polymer  preform 
matrix,  a  1-1-3  piezocomposite  with  high  hydrostatic  charge  constant  can  be  obtained.  Sherrit, 
Wiederick  and  Mukherjee  (1992)  also  described  a  1-3  type  PZT- air  piezocomposite  with 
extremely  high  effective  hydrostatic  piezoelectric  constants  and  figure  of  merit. 

4.4.1  Lateral  stress  reduction  and  the  role  of  interlayer 

The  influence  of  an  interlayer  with  a  thickness  of  0.1  mm  is  shown  in  Fig.  4.10,  where  the 
is  the  interlayer  modulus  and  T°  is  the  modulus  of  plain  Spurr  epoxy  as  given  in  Table  2.3 
(Case  1).  The  PZT-5H  rod  aspect  ratio  is  fixed  at  13.3.  The  introduction  of  a  thin  soft 
interlayer  greatly  reduces  the  lateral  stress  in  the  PZT  rod.  The  effect  of  softening  the  interlayer 
becomes  dominant  for  rod  volume  fractions  between  20%  and  45%  where  the  increases 

significantly  with  a  decrease  in  the  interlayer  stiffness. 

The  effect  of  the  lateral  stress  reduction  for  the  case  of  a  soft  interlayer  is  illustrated  in  Fig. 
4.11.  The  average  stresses  in  the  PZT  rods  are  plotted  in  Fig.  4.11  along  with  the  for 

Y(^)  =  y°  / 100  .  A  comparison  with  Fig.  4.7  reveals  that  there  is  a  large  reduction  in  the  lateral 
stress  for  the  comphant  interlayer.  Although  the  axial  stress  is  also  shghtly  decreased, 
competition  among  the  axial  stress,  the  lateral  stress  and  the  volume  fraction  of  the  PZT  rod 
gives  an  enhanced  effective  hydrostatic  piezoelectric  constant  .  This 
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Fig.  4. 10.  Variation  of  with  volume  fraction  of  PZT  rods  for  different  interlayer 
moduli. 


84 


Piezoelectric  constant,  d  (10"  mA/^) 


d  ,  Y^^^=Y®/100 

h 


- (T  +T  J/(-p) 

rr  00 

- (T  )/(-p) 

zz 


PZT  rod  volume  fraction,  v 

f 


Fig.  4. 1 1 .  Variation  of  and  the  average  stresses  in  the  rods  with  rod  volume  fraction 
for  the  case  of  a  very  soft  interlayer. 
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Normalized  average  stresses 


lateral  stress  reduction  effect  can  be  further  explained  by  plotting  the  product  of  the  average 
stress  in  the  PZT  rods,  the  rod  piezoelectric  constant  and  the  rod  volume  fraction  in  Fig.  4.12 
for  the  same  case  as  in  Fig.  4.11.  According  to  Eq.  (4.4),  the  difference  between  the  curve 
corresponding  to  the  average  axial  stress  and  the  curve  related  to  the  lateral  stress  is  exactly  the 
effective  hydrostatic  constant  d^,  which  is  also  shown  in  Fig.  4.12.  Thus,  lateral  stress 

reduction  through  the  introduction  of  a  thin,  soft  interlayer  is  also  beneficial  for  the  enhancement 
of 

4.4.2  Influence  of  the  interlayer  Poisson’s  ratio  on 
The  Poisson’s  ratio  of  the  polymer  matrix  is  also  a  very  important  parameter  in  designing  1-3 
piezocomposites  for  hydrostatic  apphcations.  If  the  Poisson's  ratio  of  the  polymer  matrix  is 
large,  the  polymer  matrix  will  be  hydrostatically  incompressible  and  large  lateral  stresses  will 
develop  in  the  piezoceramic  rods.  Previous  solutions  to  suppress  the  Poisson  effect  of  the 
polymer  matrix  rehed  on  introducing  bubbles  into  the  polymer  matrix  to  reduce  the  Poisson’s 
ratio.  Voids,  however,  will  also  reduce  the  matrix  .stiffness  and  introduce  an  undesired  bias 
pressure  dependence  to  the  hydrostatic  sensitivity  (Haun  and  Newnham,  1986;  Kim, 
Rittenmyer  and  Kahn,  1993) 

The  lateral  stress  in  the  piezoceramic  rods  can  also  be  reduced  by  controlling  the  Poisson’s 
ratio  of  a  thin  interlayer.  This  effect  is  shown  in  Fig.  4.13  for  the  case  of  =  T”  /  10.  The 
thickness  of  the  interlayer  is  fixed  at  0.1  mm  and  the  PZT-5H  rod  aspect  ratio  is  13.3.  Figure 
4.13  shows  a  comphant  interlayer  with  a  smaller  Poisson’s  ratio  can  further  reduce  the  lateral 
stress  in  the  PZT  rod  and  increases  the  of  the  1-3  piezocomposite.  The  maximum  value 

is  again  shifted  to  higher  PZT  rod  volume  fraction. 
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Interactions  between  the  axial  and  lateral  stress  versus  the  rod  volume 
fraction. 
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Fig.  4. 13.  Variation  of  with  volume  fraction  of  PZT  rods  for  different  interlayer 
Poisson’s  ratios. 


88 


4.5  Finite  Element  Simulation  of  Piezocomposites  with  Tailored  Interphase 

In  the  previous  discussion,  lateral  stress  reduction  and  axial  stress  amplification  were  not 
obtained  simultaneously.  In  order  to  improve  the  hydrostatic  performance  of  a  1-3 
piezocomposite  further,  the  interlayer  should  be  tailored  to  enhance  both  these  effects.  An 
interlayer  that  was  functionally  graded  along  the  axis  of  the  rod,  such  that  it  was  stiffer  at  the 
edges  of  the  rod  to  enhance  load  transfer  and  axial  stress  amplification  and  softer  along  the 
interior  region  of  the  rod  to  reduce  lateral  stresses,  should  be  able  to  exploit  both  effects. 
Based  on  this  concept,  a  tailored  interlayer  design  was  developed  as  shown  schematically  in 
Fig.  4.14.  Evaluation  of  the  hydrostatic  response  of  a  1-3  piezoelectric  composite  with  such  an 
interlayer  was  conducted  using  the  finite- element  method.  Because  of  the  increased  complexity 
of  this  tailored  interphase,  a  numerical  method  is  necessary  to  study  the  influence  on  the 
electromechanical  couphng  in  the  composite. 

A  composite  cyhnder  model  is  used  to  represent  a  single  piezoceramic  rod,  the  interlayer 
and  the  surrounding  matrix.  The  interlayer  consists  of  two  parts,  the  core  and  the  edge,  each 
with  distinctive  material  properties.  As  an  example,  the  interlayer  edge  was  chosen  as  the  same 
as  the  matrix.  Symmetry  in  the  axial  or  z  coordinate  is  preserved.  The  composite  model  is 
developed  with  8- node  quadrilateral  elements  (CAX8RE  elements  in  ABAQUS).  The  finite 
element  mesh  had  a  total  of  5400  elements  and  is  shown  in  Fig.  4.15.  The  length  of  the 
composite  cyhnder  model  21  is  10  mm  and  the  length  of  the  interlayer  core,  2//,  is  8  mm.  The 
interlayer  is  of  0.1  mm  thick.  The  diameter  of  the  piezoceramic  rod  is  0.75  mm.  Fig.  4.16 
shows  the  finite  mesh  near  the  end  of  the  interlayer  core.  The  finite- element  analysis  was 
conducted  using  ABAQUS  (Hibbitt  Karlsson  &  Sorenson  Inc.)  software.  The  boundary 
conditions  apphed  are  shown  schematicaUy  in  Fig.  4.17.  Material  properties  for  the 
piezoceramic  and  epoxy  matrix  are  from  Table  2.1  and  Table  2.3  (Case  1). 
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Fig.  4.14.  Schematic  of  a  1-3  piezocomposite  with  a  tailored  interlayer. 
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Fig.  4.16.  Detail  of  the  finite -element  mesh  near  the  edge  of  the  soft  interlayer. 
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The  calculated  effective  piezoelectric  constant  along  with  values  from  theoretical 
predictions  are  plotted  in  Fig.  4.18.  Finite- element  results  are  depicted  by  open  or  closed 
symbols  while  sohd  or  dashed  lines  represent  analytical  solutions.  The  composite  with  a  graded 
interlayer  has  a  broad  peak  when  compared  with  the  comphant  interlayer  case  and  has 
improved  hydrostatic  behavior  compared  with  the  no  interlayer  case.  This  broadened  peak  for 
the  effective  piezoelectric  constant  d,^  leaves  room  for  the  consideration  of  other  design 

parameters,  such  as  the  rod  volume  fraction  that  is  critical  to  the  effective  dielectric  constant  and 
the  acoustic  impedance. 

4.6  Discussion 

The  analysis  of  the  hydrostatic  performance  of  1-3  piezocomposites  demonstrated  how  to 
optimize  electromechanical  couphng  in  the  design  of  a  piezocomposite  material  for  a  specific 
application:  to  build  a  composite  stmcture  that  transfers  the  external  stresses  applied  to  the 
active  component  of  the  composite  in  a  pattern  that  most  nearly  approximates  the 
piezoceramic’s  maximally  coupled  stress  pattern.  The  micromechanics  model  developed  in 
Chapter  2  was  used  to  predict  the  effective  hydrostatic  piezoelectric  constants  d,^ ,  and  the 
figure  of  merit,  d^gf^ .  Parametric  studies  were  performed  to  assess  the  influence  of  matrix 

stiffness,  rod  aspect  ratio,  interlayer  stiffness,  and  rod  volume  fraction  on  the  hydrostatic 
piezoelectric  constants.  Increasing  the  aspect  ratio  significantly  increases  the  d^  due  to  the  axial 

stress  amplification  and  increased  load  transfer.  This  effect  saturates  for  an  aspect  ratio  of 
21  Id  >500  when  the  solution  approaches  the  values  predicted  by  a  generalized  plane  strain 
model  (infinite  aspect  ratio).  Even  though  a  large  aspect  ratio  would  increase  hydrostatic 
sensitivity,  piezocomposites  are  typically  fabricated  with  small  rod  aspect  ratios  due  to  poling 
limitations.  The  properties  of  these  piezocomposites  are  more  accurately  modeled  by  a  finite 
concentric  cylinder  model. 
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Fig.  4.18.  Prediction  of  for  a  1-3  piezocomposite  with  a  tailored  interlayer. 
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The  can  also  be  increased  by  decreasing  the  matrix  stiffness  or  introducing  a  comphant 

interlayer  around  the  PZT  rods.  Decreasing  the  matrix  stiffness  causes  axial  stress  amphfication, 
which  leads  to  higher  values.  However,  difficulties  may  exist  in  fabricating  an  entire 
composite  with  a  very  comphant  matrix.  A  comphant  interlayer  significantly  increases  the  by 

effectively  attenuating  the  lateral  stress  in  the  PZT  rods  without  decreasing  the  overah  stiffness  of 
the  composite.  A  comphant  interlayer  with  smaher  Poisson’s  ratio  can  further  increase  the  d^  of 

the  1-3  piezocomposite. 

To  improve  the  hydrostatic  performance  of  a  1-3  piezocomposite  further,  one  must  tailor 
the  interlayer  to  enhance  both  reduction  in  lateral  stress  and  axial  stress  amphfication.  An 
interlayer  that  is  functionally  graded  along  the  axis  of  the  rod,  such  that  it  is  stiffer  at  the  edges  of 
the  rod  to  enhance  load  transfer  and  axial  stress  amphfication  and  softer  along  the  interior  region 
of  the  rod  to  reduce  lateral  stresses,  could  exploit  both  effects.  Such  an  optimal  interlayer 
design  would  provide  improved  electromechanical  couphng  and  design  fiexibhity. 
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5.  SUMMARY  AND  CONCLUSIONS 


The  application  of  piezoceramic -polymer  composites  continues  to  grow  in  biomedical, 
automotive  and  other  industries.  Because  of  the  potential  of  piezocomposites  for  novel 
apphcations  it  is  important  to  improve  the  conceptual  understanding  of  these  materials,  to 
develop  practical  design  mles,  and  to  investigate  ways  of  optimizing  performance.  In  this 
dissertation  several  contributions  have  been  made  toward  these  goals. 

5.1  Summary  of  Results 

A  micromechanics  model  was  developed  for  predicting  the  local  fields  and  effective 
behavior  in  piezocomposites  with  1-3  connectivity.  Because  the  presence  of  a  thin  interlayer  or 
polymer  coating  around  the  ceramic  rods  can  influence  the  local  interaction  between  the 
piezoceramic  and  polymer  matrix  and  change  the  overall  performance  of  the  composite,  the 
model  was  designed  to  incorporate  an  interlayer  with  varying  properties.  This  composite  model 
was  also  characterized  by  its  finite  length,  chosen  to  be  equal  to  that  of  the  embedded 
piezoceramic  rods.  The  solutions  of  the  displacement,  stress,  and  electric  fields  in  the 
composite  model  were  obtained  using  two  displacement  potential  functions  and  an  electrical 
potential  function. 

The  micromechanics  model  was  used  to  predict  static  displacement  behavior  in  a  1-3 
piezocomposite.  Predictions  from  the  current  model  were  compared  with  solutions  generated 
by  simple  plane- strain  models  and  other  micromechanical  models.  The  results  demonstrated 
that  the  local  out-of-plane  displacement  of  the  piezocomposite  is  highly  dependent  on  the 
interaction  between  the  polymer  matrix,  a  thin  interlayer,  and  the  piezoelectric  ceramic.  In  order 
to  enhance  the  out-of-plane  displacement,  one  should  choose  a  passive  matrix  with  small 
Young’s  modulus.  This  result  is  similar  to  that  obtained  from  simple  models.  However,  the 
overall  stiffness  of  the  composite  is  dependent  on  the  matrix  modulus,  especially  for  lower 
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volume  fractions.  By  introducing  a  thin  compliant  interlayer,  one  can  greatly  increase  the  out-of¬ 
plane  displacement  of  the  rod  while  satisfying  the  requirement  for  overall  stiffness  of  the 
composite.  The  interlayer  essentially  serves  to  couple  or  decouple  the  rods  from  the  matrix, 
depending  on  its  stiffness. 

The  micromechanics  model  was  also  used  to  estimate  the  effective  piezoelectric  constants 
^33  and  ^33  for  the  composite.  Rod  aspect  ratio,  matrix  properties  and  interlayer  properties 
had  a  significant  effect  on  the  values  of  and  d^^.  Thus,  the  interlayer  can  be  tailored  to 
optimize  the  properties  and  increase  the  sensitivity  of  1-3  piezocomposites. 

Experimental  studies  focused  on  probing  surface  displacements  of  1-3  piezocomposites 
using  a  scanning  heterodyne  laser  interferometer.  Static  surface  displacements  of  1-3  PZT  rod- 
epoxy  samples  with  different  interphase  regions  were  measured  and  correlated  with  the  effective 
static  performance  of  the  composite.  Several  types  of  interphase  regions  were  considered. 
Coatings  with  elastic  modulus  lower  than  that  of  the  epoxy  matrix  were  apphed  to  the  rods.  The 
influence  of  a  silane  couphng  agent  was  also  investigated.  Experimental  displacement  profiles 
were  compared  with  micromechanical  predictions  using  the  finite  composite  cylinder  model. 
Micro- interferometric  measurements  as  well  as  the  theoretical  predictions  demonstrated  that  the 
out- of- plane  displacement  of  the  composite  is  highly  dependent  on  the  properties  of  the  polymer 
matrix,  a  thin  interlayer,  and  the  piezoelectric  ceramic.  The  variation  of  surface  displacement 
profile  due  to  changes  in  interphase  properties  as  predicted  using  the  composite  model  were 
compared  with  the  experimental  measurements.  The  predictions  were  in  good  agreement  with 
the  experimentally  observed  trends.  The  out-of-plane  displacement  of  the  rod  increased 
significantly  when  a  thin  comphant  interlayer  was  introduced.  A  silane  couphng  agent  improved 
the  interfacial  adhesion  and  reduced  the  nonuniformity  of  the  displacement  profile  in  the  whole 
composite  region.  Eurthermore,  the  measured  displacement  profile  can  also  be  used  to  estimate 
the  effective  piezoelectric  constants  d^^  of  the  composite. 
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The  study  of  electromechanical  coupling  in  piezocomposites  was  further  developed  by 
investigating  the  hydrostatic  performance  of  1-3  piezocomposites,  an  important  issue  in  design 
of  piezocomposites  for  low  frequency  applications.  Emphasis  vas  placed  on  determining  the 
stress  transferred  between  the  passive  matrix  and  the  active  piezoceramic  rods  and  using  this 
data  to  indicate  the  level  of  electromechanical  coupling.  The  stress  field  in  the  piezoelectric 
ceramic  under  hydrostatic  loading  was  predicted  using  the  analytical  micromechanical  model 
developed  and  a  finite  element  model  as  well.  Optimal  electromechanical  couphng  was 
achieved  when  a  certain  favorable  stress  field  was  induced  in  the  piezoceramic. 

The  key  principle  in  designing  a  piezoelectric  ceramic -polymer  composite  is  to  maximize  the 
electromechanical  coupling.  Maximum  couphng  is  achieved  by  building  a  composite  stmcture 
which  transfers  the  apphed  external  stresses  to  the  composite’s  piezoelectric  component  in  a 
pattern  that  most  nearly  approximates  the  piezoceramic’s  maximaUy  coupled  stress  patterns. 
This  investigation  has  shown  that  the  improved  hydrostatic  performance  of  1-3  piezocomposites 
comes  from  two  major  effects:  the  axial  stress  amphfication  and  the  lateral  stress  reduction  in  the 
piezoelectric  rods.  Thus,  in  designing  1-3  piezocomposites  for  better  hydrostatic  behavior,  the 
primary  goal  is  to  increase  the  axial  stress  and  to  reduce  the  lateral  stress  in  the  piezoceramic 
rods. 

The  micromechanics  model  developed  was  used  to  examine  systematicaUy  the  effective 
hydrostatic  piezoelectric  constants  and  the  figure  of  merit,  .  Parametric  studies 

were  performed  to  assess  the  influence  of  matrix  stiffness,  rod  aspect  ratio,  interlayer  stiffness, 
and  rod  volume  fraction  on  the  hydrostatic  piezoelectric  constants.  Increasing  the  aspect  ratio 
significantly  increased  the  due  to  axial  stress  amplification  and  increased  load  transfer.  This 

effect  saturates  for  2lld  >  500 ,  where  the  solution  approaches  the  values  predicted  by  a 
generahzed  plane  strain- model.  Even  though  a  large  aspect  ratio  would  increase  hydrostatic 


99 


sensitivity,  piezocomposites  are  typically  fabricated  with  small  rod  aspect  ratios  due  to  pohng 
limitations.  The  properties  of  these  piezocomposites  are  more  accurately  modeled  by  a  finite 
concentric  cyhnder  model. 

The  df^  can  also  be  increased  by  decreasing  the  matrix  stiffness  or  introducing  a  comphant 

interlayer  around  the  PZT  rods.  Decreasing  the  matrix  stiffness  causes  axial  stress  amphfication, 
which  leads  to  higher  values.  However,  difficulties  may  exist  in  fabricating  an  entire 
composite  with  a  very  comphant  matrix.  A  comphant  interlayer  significantly  increases  the  by 

effectively  attenuating  the  lateral  stress  in  the  PZT  rods  without  decreasing  the  overall  stiffness  of 
the  composite.  A  comphant  interlayer  with  smaller  Poisson’s  ratio  can  further  increase  the 

of  the  1-3  piezocomposite. 

To  improve  the  hydrostatic  performance  of  a  1-3  piezocomposite  further,  one  must  tailor 
the  interlayer  to  enhance  both  reduction  in  lateral  stress  and  axial  stress  amphfication.  It  has 
been  shown  theoretically  that  an  interlayer  mechanically  graded  along  the  axis  of  the  rod,  such 
that  it  is  stiffer  at  the  edges  of  the  rod  to  enhance  load  transfer  and  axial  stress  amphfication  and 
softer  along  the  interior  region  of  the  rod  to  reduce  lateral  stresses,  should  be  able  to  exploit 
both  effects.  This  optimal  interlayer  design  provides  improved  electromechanical  couphng  and 
design  flexibility. 

Although  the  current  work  is  focused  on  the  electromechanical  behavior  of  1-3 
piezocomposites  at  low  frequency,  the  research  results  and  conceptual  understanding  obtained 
have  importance  for  optimizing  the  design  of  a  piezocomposite  in  other  apphcations  as  well. 

5.2  Directions  for  Future  Work 

One  strength  of  this  micromechanical  analysis  is  the  abihty  to  explore  better  approaches  to 
optimization  of  the  electromechanical  behawor  of  piezocomposites.  The  method  can  be  used  to 


100 


study  the  effect  of  the  piezoceramic/polymer  interface/interphase  further.  The  results  could 
provide  a  sound  base  for  optimal  interface/interphase  design  and  control  in  piezocomposites. 

The  micromechanics  methods  can  be  applied  to  study  piezoelectric  composites  with  other 
geometrical  configurations  such  as  2-2  piezocomposites  and  composites  with  particulates.  The 
methods  can  also  be  used  to  investigate  other  coupled  phenomena  in  composite  materials,  such 
as  thermo- mechanical  and  pyroelectric  effects. 

The  micro- interferometric  method  is  ideally  suited  to  investigate  the  local  deformation  in 
piezocomposites,  especially  near  the  ceramic/polymer  interface.  Accurate  measurement  of  the 
effective  piezoelectric  constants  for  the  piezocomposites  is  stiU  a  difficulty  in  practical 
applications  of  the  composite  materials.  The  interferometer  could  be  adapted  to  measure  the 
tme  effective  piezoelectric  constant  for  a  composite. 

The  experimental  method  can  also  be  used  to  investigate  dynamic  response  of  the 
piezocomposites  at  higher  frequencies.  The  method  can  be  used  to  measure  different  types  of 
resonance  modes  of  piezocomposites  and  to  study  rod- rod  interactions.  A  significant  amount  of 
work  remains  to  be  done  in  this  direction. 
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APPENDIX 


The  coefficients  in  Eq.  (2.30)  are  given  by 
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The  functions  where  7=1, 2,. ..24,  in  Eqs.  (2.40)  to  (2.47)  can  be  expressed  in  terms  of 
modified  Bessel  functions  of  the  first  kind  /g,  /j  and  elastic  constants  of  the  piezoelectric 
ceramic  as  follows: 
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In  region  2  and  region  3  the  complementary  solution  for  the  displacements  and  stresses  are 
of  the  form 
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where  F^'^  is  Love’s  stress  function  and  the  superscript  i=2,3  for  the  interlayer  and  matrix 
region  respectively.  The  functions  and  where  i=2,3  and 7=1,2,...  12,  in  Eqs  (2.52)  to 

(2.58)  can  be  expressed  in  terms  of  modified  Bessel  functions  of  the  first  and  second  kind  Iq, 
/j ,  Kq,  Fj  and  the  elastic  constants  as  follows: 
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